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The present work includes some of the author's 
original researches on the integer solutions of equations 
and linear systems: 

1. The notion of "general integer solution" of a 
linear equation with two unknowns is extended to linear 




2 



equations with n unknowns and then, to linear systems. 

2. The proprieties of the general integer solution 
are determined (both of a linear and system) . 

3. Seven original integer algorithms (two for 
linear equations and five for linear systems) are 
presented. The algorithms are strictly demonstrated and 
an example for each of them is given. These algorithms 
can be easily introduced in a computer. 



INTEGER SOLUTIONS OF LINEAR EQUATIONS 



Definitions and properties of the integer solution of 
linear equations (l.e) : 

Let l.e: 

n 

(1) £ aiXi = b, with not all a± null and b from Z. 

i=l 

Again, let h e N, and the functions f± : Z h -> Z, 
i = 1 , n . 

Definition 1 

o 

xi = xi, i = l,n, is the particular integer solution 

o n o 

of equation (1) , if all xi e Z and £ aiXi = b. 

i=l 
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Definition 2 

xi = fi(ki, ... kh) , i = l,n, is the general integer 
solution of equation (1) if: 



n 

(a) 2 aj.fi (ki, ... k h ) = b V (ki, ... k h ) e Z h , 

i=l 



0 0 

(b) Irrespectively of f± (xi, ... x n ) 

o 

particular integer solution for (1) (ki, . 

ooo 

that Xi = fi (ki, ... k h ) for all i = l,n. 

We will further see that the general integer solution 
can be expressed by linear functions. 

h 

We consider for 1 < i < n, the functions fi = 2 

j = l 

Cijkj + di with all cij, di e Z. 



there is a 



o 

k h ) 



Z h so 



Definition 3 

A = ( Cij ) i, j is referred to as the matrix associated 
to the general solution of equation (1) . 

Definition 4 

The integers ki, . . . k s , 1 < s < h, are independent 
if all the corresponding column vectors of matrix A are 
linearly independent. 
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Definition 5 

An integer solution is s - times undetermined if the 
maximal number of independent parameters is s. 

Theorem 1 . The general integer solution of equation 
(1) is (n-1) - times undetermined. 

Proof 

We suppose that the particular integer solution is of 
the form: 

r 

(2) xi = X iiePe + Vi, i = l,n, with all u ie , VieZ, 
e=l 

P e = are parameters from Z, while a < r < n - 1. 

o o 

Let (xi, . . . x n ) a general integer solution of equation 
(1) (we are not interested in the case when the equation 
does not accept integer solution) . The solution 

o 

Xj = a n lj + Xj for j = l,n - 1 
n-1 o 

Xn = - ( X ajkj - Xn) 

j=l 

is (n-1) - times undetermined (it can be easily checked 
that the order of the associated matrix is n - 1) . Hence, 
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there are n - 1 undetermined solutions. Let, in the 
general case, a solution n - 1 times be undetermined: 



n-1 

Xi = £ Cijkj + di i = l,n, with all cij, d± e Z . 

j=l 



The case when b = 0 . 



n n 



Then £ a±Xi = 0 . 


It follows £ a±xi 


= 


i=l 


i=l 





n 


n-1 




n 


n-1 


n 


= £ a± 


(£ Cijkj + di ) 


= 


£ai 


£ Cijkj + 


£ aidi 


i = l 


j=l 




i=l 


j=l 


i=l 



n 

For kj =1 and kj = 0, j ... jo, it follows £ aiCij = 0. 
o i = 1 o 

Let the homogenous linear system of n equations with n 
unknowns be: 

n 

£ XiCij = 0, j = 1 , n - 1; 
i=l 

n 

£ xidi = 0 

i=l 

which, obviously accepts the solution xi = a±, i = l,n. 



different from the trivial one. Hence the determinant of 
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the system is zero, i.e., the vectors Cj = (Cij, ... Cnj) 11 , 

j = l,n - 1, D = (di, ..., d n ) 1 are linearly dependent. 

But the solution being n - 1 times undetermined it follows 

that Cj, j = l,n - 1, are linearly independent. Then (Ci, 

. . . C n -i) determines a free sub-module Z of the order n - 1 

in Z n of solutions for the given equation. 

Let us see what can be obtained from (2) . We have: 

n n r 

0 = £ aiXi = £ a± ( £ Ui e P e + vi) . As above, 

i=l i=l e=l 



we obtain: 

n r 

£ aiVi = 0 and £ aiUi e = 0; and similarly, the 
i=l e=l 0 

vectors Uh = (uih, . .., u n h) are linearly independent, 

h = l,r. Uh, h = l,r, and V = (v z , ..., v n ) are 
particular integer solutions of the homogenous linear 
equation . 

Subcase (al ) 

Uh, h = l,r, are linearly dependent. It gives 
{Ui, ..., U r } = the free submodule of order r in z n of 
solutions of the equation. Hence, there are solutions 
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from {Vi, V n -i } which are not from {Ui, U r }, this 

contradicts the fact that (2) is the general integer 
solution . 



Subcase (a2 ) 

Uh, h = l,r, V are linearly independent. Then, 

{Ui, U r } + V is a linear variety of dimension < n - 1 

= dim {Vi, Vn-i } and the conclusion can be similarly 

drawn. The case when b 0. 

n n n-1 

So, £ aixi = b. Then £ ai ( £ Cijkj + di) = 

i=l i=l j=l 

n-1 n n 

= £ (£ aiCij) kj + £ aidi = b, V (ki, ..., k n -i, 

j=l i=l i=l 

n 

0)e Z n_1 . As in the previous case we get £ aidi = b and 

i=l 

n 

£ aiCij = 0, Vj = l,n - 1. The vectors Cj = 
i=l 



= (cij, ..., Cnj) 1 , j = l,n - 1, are linearly independent 
because the solution is n - 1 times undetermined. 
Conversely, if Ci, ... C n -i, D (where D = (di, ..., dn) 1 ) 

n-1 

were linearly dependent, it would mean that D = £ SjCj, 

j = l 




n 

with all Sj scalar; it would also mean that b = £ aidi 

i=l 

n n-1 n-1 n 

= £ a± ( £ SjCij) = £ Sj ( £ aiCij) = 0. Impossibly, 

i=l j=l j=l i=l 

(3) then {Ci, C n -i} + D is a linear variety. 

Let us see what we can obtain from (2) . We have: 

n nr r n 

b = £ aixi = £ a± ( £ u ie P e + vi) = £ ( £aiU ie )Pe 

i=l i=l e=l e=l i=l 

n 

+ £ aivi 

i=l 

n n 

and, similarly: £aivi = b and £ aiUi e =0, V e = l,r, 
i=l i=l 

respectively. The vectors U e = (ui e , . .., u ne ) t /- e = 

l,r, are linearly independent because the solution is r 
- times undetermined. 

A procedure like that applied in (3) gives that Ui, 
..., U r , V are linearly independent, where V = (vi, ..., 
Vn) 1 . Then {Ui, ... U r } + V = a linear variety = free 
sub-module of order r < n - 1. That is, we can find 
vectors from Ci, ..., C n -i + D which are not from {Ui, 



• • • r 



U r } + V, contradicting the "general" characteristic 
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of the integer number solution. Hence, the general 
integer solution is n - 1 times undetermined. 

Theorem 2 . The general integer solution of the 

n 

homogenous linear equation £ aiXi = 0 (all ai e Z \ {0}) 

i=l 

can be written under the form: 



n-1 

(4) xi = £ Cijkj, i = l,n, (with di = ... = d n = 0) . 

j=l 

Definition 6 . This is called the standard form of 
the general integer solution of a homogenous linear 
equation . 



Proof 



We consider the general integer solution under the 
n-1 

form: xi = £ c±jPj + di, i = l,n, with not all di = 

j=l 

0. We show that it can be written under the form (4) . 
The homogenous equation admits the trivial solution 

xi = 0, for 

o o n-1 o 

i = l,n. There is (pi, ..., p n -i) e Z n_1 so that £ CijPj + 

j = l 

+ di =0, V i = l,n. Substituting: Pj = kj + Pj, 



j = l,n - 1, in the form from the beginning of the 
proof we obtain form (4) . We have to 



mention that the substitution does not diminish the 
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degree of 

generality as Pj e Z because kj e Z, j e l,n - 1. 

Theorem 3 . The general integer solution of a non- 
homogeneous linear equation of its associated homogenous 
linear equation + any particular integer solution of the 
non-homogeneous linear equation. 



Proof 

n-1 

Let xi = £ Cijkj, i = l,n, be the general integer 

j=l 

solution of the associated homogenous linear equation 
and. 



again, let xi = Vi, i = l,n, be a particular integer 
solution 

of the non-homogeneous linear equation. Then, 
n-1 

Xi £ Cijkj + 

j=l 

+ Vi, i = l,n, is the general integer solution of the 

n 

non-homogeneous linear equation. Actually, £ aiXi = 

i=l 

n n-1 n n-1 n 

= £ a ± ( £ Cijkj + Vi) = £ ai ( £ Cijkj) + £ aiVi 

i = l j = l i=l j = l i=l 

= b; 

o 

if Xi = Xi, i = l,n, is a particular integer solution of 
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the non-homogeneous linear equation, then xi = x ± - vi, 

i = l,n, is a particular integer solution of the 
homogenous linear 

o o 

equation; hence, there is (kg, ..., k n -i) e Z n_1 so that 

n— 1 o o n— 1 o o 

E Cijkj = xi - vi, Vi = l,n, i.e., E Cijkj + Vi = Xi,Vi = 
j=l j=l 

l,n, what had to be proven. 

n-1 

Theorem 4. If x x = E Ci-jk-j, i = l,n, is the 

j=1 

general integer solution of a homogenous linear equation 

(Cij, . . . , Cn j ) 1, 

V j = 1, n -1 . 

The demonstration is made by reductio ad absurdum. 
If 3 jo, 1 < j 0 < n - 1, 

so that (cii , ..., Cni ) = dj ± +1, 

o o o — 



then Cii 

o 



Cij dij 
0 0 



with (Cii 

o 



c n j ) =1, Vi =1, n. 

o 



But Xi = 
solution 



'ID, 



as 



, i = l,n, represents a particular 



n 




n 


n 


E aiXi 


= 


E aiCij = 1/dj ! 


E 


i=l 




i=l 


i=l 



integer 

aiCii = 0 

o 



(because Xi = Cij , i = l,n, is a particular integer 

solution from the general integer solution by 
introducing kj =1 and 
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o 



kj = 0, j jo) . But the particular integer solution 



xi = Ci-j , i = l,n, cannot be obtained, by introducing 

o 

whole number parameters (as it should) , from the general 
integer solution, as, from the linear system of n 
equations and n - 1 unknowns, which is compatible, we 
obtain : 



n ■ . 

Xi S Cijkj Cij dj kj Cij , i 1, n, 

j _ ooo o 

j*jo 

Leaving aside the last equat ion--which is a linear 
combination of the other n - 1 equations, a Kramerian 
system is obtained. It follows up: 



I Cn ... Cij ... Ci,n-1 ! 

I ° I 

I Cn-1,1 • • • Cn-li •••Cn-l,n-l i 

o 1 

kj = = e Z . 

° ’ dj o 

I Cn ... Cij dj ... Ci,n-1 I 

! ■ 0 0 I 

I Cn-l, 1 • Cn-1, j dj . Cn-l,n-l I 
0 0 

The assumption is false and thus the proof results. 

Theorem 5 . Considering the equation (1) with 

(ai, ..., a n ) =1, b = 0, and the general integer 
solution 

n-1 

x i ^ Cijkj, i 1 , n, then ( a i , ..., a i— i , a i+i , ..., 

j = l 

a n ) - (Cil, ..., Cin-i) , Vi 0 l,n. The demonstration is 
made by 




double divisibility, 
but 



Let i 



1 < i 0 < n, be arbitrary 



0 r 



fixed. xi 



n-1 

£ Ci jkj . Consider the equation £ 



= -ai xi . We have shown that xi = c±i, i = l,n is a 

oo 

particular integer solution irrespective of j, a < j 
< n - 1. The equation 6 aixi = -a± Ci j, obviously, 



admits the integer solution xi = Cij, i ... io- Then 

(a i , • • • , ai— i, ai +i, . . . , an) divrdes ai Ci j , but, 

we 

oo oo 

have assumed that (ai, ..., a n ) - 1, it follows that 

(ai, . . . , ai -i, ai +i, . . . , a n ) I Ci i, irrespective of j . 

0 0 0 

Hence (ai, . . . , ai — i, ai — i, . . . , an) I (Ci i, . . . , Ci n-i) , 

oo oo 



Vi 0 l,n, and the divisibility in one sense was 
proven . 



Inverse Divisibility 

Let us suppose the contrary and say that > ii 0 



0 l,n for which (ai, ..., ai -i, ai +i, ..., a n ) - di i ... 

2 ~ 

ii i 

- (ci i, . • ., Ci n -i ) ; we have considered di i and di 2 

11 11 



aix 



< 



as 



. di 



without restricting the generality. di 1 | di 2 according 
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to 

i i 

the first part of the demonstration. Hence, > d 0 Z so 
that 

d± 2 = d ! d± i, | d | ... 1 . 

i i 



n-1 

Xi = 6 Ci jkj = d ! di 

aiXi = 



= -an 



Xi 



0 

i-ii 



aiXi 



-a. 



n-1 n 

6 Ci jkj,* 6 aiXi = 0 Y 6 

j = l 1 i=l i-ii 



n-1 

d ! di i 6 Ci jkj, where 



( Ci i , 
1 



• • • , Ci n-1 ) 1 

1 



The nonhomogeneous linear equation 6 aiXi = -ai di 
admits i...ii 

integer solution because ai di is divisible by 

o 

( a. i , . . . , a. i — i , a. i+i, . . . , a. n ) • li "t i i , i ... i i ]o 
its 

i i 

particular integer solution. It follows that the 
equation 

n o 

6 aiXi = 0 admits the particular solution Xi = Xi, 
i=l 

i ... ii, Xi = di , which is written as (5) . We show that 
(5) 



i 



i 
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cannot be obtained from the general solution by integer 
number parameters : 



n-1 o 

6 Cijkj = xi, i ... ii; 

j = l 

n-1 

d ! di i 6 Ciiki = di i (6) 

1 j=l 



But equation (6) does not admit integer solution because 
d ! di i di i thus, contradicting, thus, the "general" 

characteristic of the integer solution. 

As a conclusion we can write: 

Theorem 6. Let the homogenous linear equation be: 



n 

6 aiXi = 0 with all a ± 0 Z \ {0}, and (a lf . .., a n ) - 

1 . 

i = l 

h 

Let Xi = 6 Cijkj, i = l,n be, with all c±j 0 Z, all kj 

j=l 

whole integer parameters and h 0 N, a general integer 
solution of the equation. Then, 

1° the solution is n - 1 times undetermined; 



2° Vj = l,n - 1 we have (cij, ..., c n j) - 1; 

3° Vi = l,n we have (cu, ..., Ci n -i) - (ai, ..., ai-i. 
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a i+ i, a n ) . The proof results from Theorems 1, 4 and 

5 . 

Note 1 . The only equation of the form (1) which is 
n - times undetermined is the trivial equation 0 ! xi + 

+ . . . + 0 ! x n = 0 . 

Note 2 . The converse of theorem 6 is not true. 
Counterexample : 



Xi = 


- ki 


+ k 2 




x 2 = 


5kj 


+ 3k 2 




X 3 = 


7k x 


- k 2 , ki, k 2 0 Z 


(7 



is not the general integer solution of the equation 



-13xi + 3 x 2 



4x 3 = 0 



although the solution (7) verifies the points 1°, 2° and 



of theorem 6. (1, 7, 2) is the particular integer 

solution 

of (8) but cannot be obtained by introducing integer 
number 

parameters in (7) because from 



-kg + k 2 — 1 



5ki + 3k 2 



7 
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7ki - k 2 =2 



follows that k = 1/2 0 Z and k = 3/2 0 Z (unique roots) . 



Reference : 

[1] Smarandache, Florentin, “Whole number solutions of 
linear equations and systems”, M. Sc. Thesis, 1979, 
University of Craiova (under the supervision of 
Assoc. Prof. Dr. Alexandru Dinca) . 
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AN INTEGER NUMBER ALGORITHM 
TO SOLVE LINEAR EQUATIONS 

An algorithm is given which ascertains whether a 
linear equation admits integer number solutions or not; if 
it does, the general integer solution is determined. 

Input 

A linear equation aiXi + ... + a n x n = b with ai b 0 Z, 

xi being integer number unknowns, i = l,n and not all 
ai = 0 . 

Output 

Decision on the integer solution of this equation; if 
the equation has solutions in Z, its general solution is 
obtained . 

Method 

Step 1. Calculate d = (ai, . . ., a n ) . 

Step 2. If d/b then "the equation has integer 
solution"; go on to Step 3. If d/b then "the equation 
does not admit integer solution"; stop. 



Step 3. Consider h: 



1. If | d | ... 1 divides the 
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equation by d; consider ai : = a±/d, i = l,n, b: = b/d. 



Step 4 . 


Calculate 


a = min | 

a s ...O 


a s | and determine an i so 


that a± = a . 










Step 5 . 


If a ... 


1, 


go on to 


Step 7 . 


Step 6. 


If a = 


If 


then : 





(A) xi = -(aixi + ... + ai-ixi-i + ai+iXi+i + 
+ ... a n x n - b) ! ai 



(B) Substitute the value of Xi in the values of 
the other determined unknowns . 

(C) Substitute integer number parameters for 

all the variables of the unknown values in 
the right term: ki, k 2 , . .., k n -2 and k n -i, 

respectively . 

(D) Write down the general solution thus 
determined; stop. 

Step 7. Write down all aj, j ... i and b under the 

f orm : 

aj = aiqj + rj 

aj b 

b = a±q + r, where qj = -- , g = 

9. j_ 3. y 

Step 8 . Write x ± = - qiXi - ... - q±-iXi-i - qi+iXi+i - 
- . . . - q n x n + q - th. Substitute the value of xi in the 
values of the other determined unknowns. 
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Step 9. ai: = 4i 
Consider 

ai-1 : = r 

a.i+i • = ri+i 



and ai : = -ai 
b : = r 
Xil = t h 

h : = h + 1 



a n : = r n and go back to Step 4. 

Lemma 1 . The previous algorithm is finite. 



Proof 

Let the initial linear equation be aiXi + ... + a n x n = 
= b with not all a± = 0; it is considered that min | a s I = 

a s ...O 

ai ... 1 (if not, it is renumbered) . Following the 
algorithm, once we pass from this initial equation to a 
new equation: aiti + a 2 X 2 + ... + a n x n = b', with | ai | < 

< | ai | for i = 2,n, |b' | < |b| and ai = - ai. It 

follows that min | a s I < min | a s | . We continue similarly 

a'...0 a s ...l 

S 

and after a finite number of steps we get, at Step 4, at 
a: = 1 (as, every time, at this step the actual a is 
strictly smaller than the previous a, according to the 
former note) and in this case the algorithm terminates. 
Lemma 2 . Let the linear equation be: (25) aiXi + 

+ a 2 x 2 + . . . + a n x n = b with min | a s I = ai and the 

a s ...O 
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equation: (26) - aiti + r 2 x 2 + . . . + r n x n = r with ti = 



= - xi - q 2 x 2 - ... - q n x n + q where r ; = a ; - aiqi, i = 

ai b 

= 2,n, r = b - aiq while q± = -- , r = -- . Then 

3 . 3.1 

0 0 0 

xi = xi, x 2 = x 2 , . . . , x„ = x n is particular solution of 

0 0 

equation (25) if and only if ti = ti = - xi - q 2 x 2 - . . . - 
o o 

q n x n + q, x 2 , . . . , x n = x n is a particular solution of 

equation (26). 



Proof 

oo o 

xi = xi, x 2 = x 2 , ..., x n = x n is a particular solution 

0 0 0 0 
of equation (25) ] aiXi + a 2 x 2 + ... + a n x n = b ] aiXi + 











0 




0 


0 


+ 


(r 2 


+ 


a 2 q 2 ) 


X 2 + ... + (r n 


+ aiq n 


) x n = aiq + r ] 


r 2 x 2 + 








0 


0 0 




0 


0 


+ 


. . . 


+ 


r n X n 


- ai (-Xi -q 2 x 2 


- ... 


- q n x n + q) = r 


] - aiti + 




0 






0 


0 


0 


0 


+ 


r 2 x 2 


+ 


. . . 


+ r n x n = r ] ti 


= ti. 


X 2 X 2 f • • . r X n 


= x n is a 



particular solution of equation (26) . 

Lemma 3 . Xi = cuki + . . . + Ci n -ik n -i + d±, i = l,n is 
the general solution of equation (25) if and only if: 
(28) 



ti — - (Cn + q 2 c 2 i + ... + q n c n i)ki - ... - (Ci n -i + 

+ q 2 C 2n -i + . . . + q n c nn-1 )Kn - (di + q 2 d 2 + ... + q n dn) + 

+ q, Xj = Ccjiki + ... + C;n--ik n + dj, j = 2 , n 
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is a general solution for equation (26) . 

Proof 

ooo o o 

ti = ti = - xi - q 2 x 2 - ... ~ q n x n + q, x 2 = x 2 , . . 

o 

. .., x n = x n is a particular solution of the equation (25) 

0 0 0 

] xi = xi, x 2 = x 2 , . . . , x n = x n is a particular solution of 

o o 

equation (26) ] > kq = ki 0 Z, k n = k n 0 Z so that 

o oo 

xi = ciiki + ... + Cin-ikn-i + di = x ±/ i = 1, n ] > ki = 

oo oo 

= ki > Z, . . ., k n k n 0 Z so that xi = Cuki + . . . + Ci n -.ik n -i. + 

o o 

+ di = Xi, i = 2 , n and tj = - (c u + q 2 c 2i + ... + q n c n i)ki - 

o 

- ... - (Cin-i + q 2 c 2n -i + ... + q n Cnn-i) k n -i - (di + q 2 d 2 + 

oo oo 

+ . . . + q n d n ) + q = - Xi - q 2 x 2 - ... - q n x n + q = ti. 

Lemma 4 . The linear equation (29) aiXi + a 2 x 2 + ... + 
+ a n x n = b with | ai | =1 has the general solution: (30) 

Xi = - (a 2 k 2 + ... + a n k n - b) a lr 

Xi = ki 0 Z, 

i = 2 , n . 



Proof 

oo o 

Let xi = Xi, x 2 = x 2 , . . ., x n = x n be a particular 

o o 

solution of the equation (29). > k 2 = x 2 , ..., k n = k n so 



o 



o 



o 



o 
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that xi = - (a 2 x 2 + . . . + a n x n - b) ai = xi, x 2 = x 2 , . . . , 

o 

• • • r X n X n . 

Lemma 5 . Let the linear equation be aiXi + a 2 x 2 + 

+ ... + a n x n = b, with min | a s | = ai and ai = aiqi, i = 2,n. 

a s -0 

Then, the general solution of the equation is: 

xi = - (q 2 k 2 + . . . + q n k n - q) , 

Xi = ki 0 Z, 

i = 2 , n . 



Proof 

Dividing the equation by ai the conditions of Lemma 4 
are met . 

Theorem of Correctness . The preceding algorithm 
correctly calculates the general solution of the linear 
equation aiXi + . . . + a n x n = b with not all a± = 0 . 

Proof 

The algorithm is finite according to Lemma 1. The 
correctness of steps 1, 2, and 3 is obvious. At step 4 

there is always min | a E | as not all ai = 0 . The 

a s ...O 



correctness of substep 6 (A) results from Lemmas 4 and 5 
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respectively. This algorithm represents a procedure of 
obtaining the general solution of the initial equation by 
means of the general solutions of the linear equation 
obtained after the algorithm was followed several times 
(according to Lemmas 2 and 3) ; from Lemma 3 it follows 
that to obtain the general solution of an initial linear 
equation is equivalent to calculate the general solution 
of an equation at step 6 (A) , equations whose general 
solution is given in algorithm (according to Lemmas 4 and 
5) . The theorem of correctness has been fully proven. 

Note . At step 4 of the algorithm we consider 

a: = min | a s I so that the number of iterations be as small 

a s ...O 

as possible. The algorithm works if we consider 

a : = | ai | ... m ax | a s I but it takes longer. The algorithm 

s=l , n 

can be introduced in the computer. 

Application 

Calculate the integer solution of the equation: 

6xi - 12x 2 - 8x 3 + 22x 4 = 14 

Solution 

The former algorhythmus is applied. 

1 . 



(6 



12 



8 , 22 ) 



2 
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2. 2 | 14 so that the solution of the equation is in 
Z 

3. h: = 1 

|2| ... 1; dividing the equation by 2 we get: 





3xi 


= 6x 2 


- 4x 3 + llx 4 = 7 


4 . 


a : 


= min 


{ 1 3 | 


, |-6| , |-4|, |11|} = 


5 . 


a .. 


. 1 






7 . 


-6 


= 3. ( 


+ 

CM 

1 


0 




-4 


= 3. ( 


+ 

CM 

1 


2 




11 


= 3.3 


+ 2 






7 


= 3.2 


+ 1 




8 . 


Xi ; 


= 2x 2 


+ 2x 3 


- 3x 4 + 2 - ti 


9. 


a 2 : 


= 0 


ai : 


= -3 




a 3 : 


= 2 


b : 


= 1 




a4 : 


= 2 


xi : 


= ti 








h: 


= 2 


4 . 


We 


have 


a new 


equation : 




-3t 


i + 0 


! x 2 


+ 2x 3 + 2x 4 = 1, 




a = 


: min 


{ 1-3 1 


, |2 | , |2 | }= 2, and 




i = 


= 3 






5 . 


a .. 


. 1 






7 . 


-3 


= 2 ! 


(-2) 


+ 1 



0 



2 ! 0 + 0 
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2 


= 2 ! 


1 + 


0 




1 


= 2 ! 


0 + 


0 


8 . 

value 


x 3 

of 


= 2ti 
x 3 in 


- 0 
the 


! x 2 - x 4 + 0 - t 2 . Substituting 
value determined for x 4 we get: 




Xi 


= 2x 2 


- 5x 


4 + 3ti — 2t 2 + 2 


9. 


ai : 


= 1 


a 3 : 


= -2 




a 2 : 


= 0 


b : 


= 1 




a 4 : 


= 0 


x 3 : 


= t 2 








h: 


= 3 


4 . 


We 


have 


obtained the equation: 




1 ! 


t 2 + 


0x 2 


- 2t 2 + 0 ! x 4 = 1, 



a = 1, and 
i = 1 



6. (A) ti = - (0 ! x 2 - 2t 2 + 0 ! x 4 - 1) ! 1 = 

2 1 2 + 1 

(B) Substituting the value of ti in the values 
of xi and x 3 previously determined, we get: 
xi = 2 x 2 - 5x 4 + 4t 2 + 5 and 

x 3 = — x 4 + 3t 2 + 2 

(C) x 2 : = ki, x 4 : = k 2 , t 2 = k 3 , k 4 , k 2 , k 3 0 Z 

(D) The general solution of the initial 
equation is: 

Xi = 2ki - 5k 2 + 4k 3 + 5 
x 2 = ki 

x 3 = -k 2 + 3k 3 + 2 
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x 4 = k 2 

ki, k 2 , k 3 are parameters 0 Z 

Reference 

[1] Smarandache, Florentin, Whole number solutions of 
equations and systems of equat ions--diploma paper. 
University of Craiova, 1979. 
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ANOTHER INTEGER NUMBER ALGORITHM TO SOLVE 
LINEAR EQUATIONS (USING CONGRUENCY) 

In the present part a new integer number algorithm 
for linear equations is presented. This is more "rapid" 
than W. Sierpinski's presented in [1] in the sense that it 
reaches the general solution after a smaller number of 
iterations. Its correctness will be strictly 
demonstrated . 



INTEGER NUMBER ALGORITHM TO SOLVE 
LINEAR EQUATIONS 

Let us consider the equation (1); 

(the case ai, b are in Q, 

i = l,n, is reduced to the case (1) by bringing to the 
same denominator and eliminating the denominators) . 

Let d = (ai, ..., a n ) . If d|b, then the equation does not 
admit integer solutions, while if d|b, the equation admits 
integer solutions (according to a well-known theorem from 
the theory of numbers) . 

If the equation accepts solutions and d ... 1, we 
divide the equation by d. Then, we can agree that d = 1 
(we do not make any restriction if we consider the maximal 
co-divisor positive) . 
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(a) Also, if all a± = 0, the equation is trivial; it 
admits the general integer solution x± = kg 0 Z, i = l,n 
when b = 0 (the only case when the general integer 
solution is n - times undetermined) and does not have 
solutions when b ... 0. 

(b) If > i, 1 < i < n, so that a± = + 1, then the 
general integer solution is: 

n 

xi = -a± ( 6 ajkj - b) and x s = k s 0 Z, s 0 

j=l 

j-i 

0 {1, . . n } \ { i } 

The proof of this assertion was give in [4] . All these 
cases being trivial, we will leave them aside. The 
following algorithm can be written: 

Input 

A linear equation: (2) 

n 

6 aixi = b, ai, b Z, ai ... + 1, i = l,n 
i=l 

with not all a± = 0 and (ai, ..., a n ) = 1. 
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Output 

The general solution of the equation. 

Method 

1 . h : = 1 , p : = 1 

2. Calculate min { | r | , r / ai (mod aj) , 

l<i, j<n 

r ... 0, | r | < | aj | } 

and determine r and the paid (i,j) for which 
this minimum can be obtained (when there are 
more possibilities we have to choose one of 
them) . 

3. If | r | ... 1 go on to step 4. 

If | r | = 1 , then 

n 

xi: = r(-ajt h - 6 a s x s + b) 

s=l 

so { i , j } 



x j : 


ai-r 

= r (a t + 


n r-ai 

! 6 a s x s + 




aj 


s=l aj 

so { i , j } 



(A) Substitute the values thus determined of 



these unknowns in all the relations (p) 
p = 1, 2 ... (if possible) . 
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the 



(B) From the last relation (p) obtained in 

algorhythmus substitute in all relations: 



(p - 1) , (p - 2) , . . . , (1) . 



(C) 



parameters 

with 



Every relation, starting in order from 

(p - 1) should be applied the same 

procedure as in (B) ; then (p - 2), 

(3), respectively. 

Write the values of the unknowns xi, 

i = l,n from the initial equation writing 
the corresponding integer number 

from the right term of these unknowns 



ki, . . . , k n -i) , STOP . 
4. Write the relation (p) : 



ai-r 

Xj = t h - xi; 

aj 

5. Consider Xj : = t h h: = h + 1 

a±:=r p : = p + 1 

The other coefficients and variables remaining 
unchanged and go back to step 2 . 
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The Correctness of the Algorithm 

Let us consider linear equation (2) . Under these 

conditions, the following proprieties exist: 

Lemma 1 . The set M = { | r | , r / ai (mod a-j) , 

0 < | r | < | aj | } has a minimum. 

Proof 

Obviously M d N* and M is finite because the equation 
has a finite number of coefficients: n and considering 

all the possible combinations of these, by twos, there is 

2 

the maximum AR n (arranged with repetition) = n elements. 
Let us show, by reductio ad absurdum, that M ... , 

M = ] a± = 0 (mod aj) Vi, j 0 l,n. Hence aj = 

= 0 (mod ai) , Vi, j 0 l,n. Or this is possible only 

when | a ± | = | a j | , Vi, j 0 l,n, which is equivalent to 

(ai, ..., a n ) = ai, Vi 0 l,n. But (ai, ..., a n ) = 1 
according to a restriction from the assumption. It 

follows that | ai | = l,n, Vi 0 l,n a fact which contradicts 
the other restrictions of the assumption. 



M ... 0 and finite, it follows that M has a minimum. 
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Lemma 2 . If |r| = min M then |r| < | a± | , Vi 0 

l<i, j<n 



0 1 , n . 



Proof 

We assume, conversely, that >i 0 , 1 < i 0 < n, so that 

I r | > | ai |. Then |r| > min { | a d | } = | aj | ... 1, 1 < jo < 

° l<j<n ° 

< n. Let a p , 1 < p 0 < n so that | a p | > | aj | and a p is 

— o — — o o o 

o 

not divided by aj . There is such a coefficient in the 

equation as | aj | is the minimum and not all the 

coefficients are equal among themselves (conversely, it 

would mean that (ai, ..., a n ) = ai = + 1, which is against 

the hypothesis) and, again, of the coefficients whose 

module is greater than | aij | not all can be divided by aj 

(conversely, it would similarly result that (ai, ..., 
a n ) = 



= aj ... + 1) . We write 1 
o — 


a p 

0 


/aj 

0 


] = qo 


0 Z (the whole 


number part) , and r = 


a p 

0 


- 


qoaj 


0 Z. 

0 


We have a p / r 

0 


(mod aj ) and 0 < | r 0 1 

0 


< 


1 aj 


1 < 

0 


1 ai I 

0 


< | r | . Thus, we 


have found a r 0 with | 


r 0 1 


< 


1 r | 


which 


contradicts the 


definition of minimum 


given 


to 


| r | . 


Contrary to the 



assumption. Thus, |r| < | a± | , V i 0 l,n. 

Lemma 3 . If |r| = min M = 1, for the pair of indices 
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( i , j ) , then : 



n 



Xi = r 


( a jth 


0 


a s k s 


+ b) 






S = 1 

so { i. 


j} 








ai-r 


n 


r-ai 


x j = r 


(ait h + 


— 


6 


a s k s + b) 






aj 


S = 1 

so { i , j } 


aj 


CO 

Ai 

II 

CO 

X 


0 Z, s 


0 {1, . 


. . . , n} 


\ { i /• j } 



is the general integer solution of equation (2) . 



Proof 

o 

Let x e = x e , e = l,n be a particular integer solution 

o 

of equation (2). Then > k s = x s 0 Z, s 0 {1, n} \ 

o ai-r 0 

\ { i , j } and t h = Xj + xi 0 Z (because ai-r = ;aj) , so 

a j 

that : 



o ai-r 

Xi = r -aj (Xj + Xi 

ai 



n 0 

6 



a s x s + b 
s=l 

so { i , j } 



= Xi 



o ai-r o 

Xj = r -aj (Xj + Xi) 

ai 



ai-r n 0 






0 a s x s + 

s=l 
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so { i , j } 



r-ai o 

+ b = Xj 

a j 



0 

and x s = k s = x s , s 0 {1, n} \ {i,j}. 



Lemma 4 . Let |r| ... 1 and (i,j) be the pair of 
indices for which this minimum can be obtained. Again, 
let the system of linear equation be: 



n 

a-jth + rxi + 6 a s x s = b 

s=l 

s...{ i, j } 

(3) 

a±-r 

th = Xj + Xi 

aj 



o 

Then, x e = x e , e = l,n is a particular integer solution 

o 

for (2) if and only if x e = x e , e 0 {1, ..., n} \ {j} and 

o o ai-r 

th = th = Xj + xi is the particular integer solution 

a j 

of (3) . 



Proof 

x e = x e , e = l,n is a particular integer solution for 

( 2 ) 
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no __ n oo ai-r 0 

] 6 a e x e = b ] 6 a s x s + aj (xj + xi) + 

e=l s=l aj 

s...{ i, j } 

o 

+ rxi = b 



o o n ^ o oo 

] ajth + rxi + 6 a s x s = b and th = Xj + 

s=l 

s-.{ i, j } 



a±-r o 

+ xi 0 Z 

aj 

0 

] x e = x e , e 0 {1,2, ..., n} \ {j} 

0 

and th = th is a particular integer solution for (3) . 

Lemma 5 . The former algorithm is finite. 

Proof 

When | r | =1 the algorithm stops at step 3. We will 
discuss the case when |r| ...1. According to the 
definition of r, r| ON*. We show that the row of r - s 
successively obtained by following the algorithm several 
times is strictly decreasing to 1. Let ri be the first 
obtained by following the algorithm one time. | ri | ... 1, 



go on to steps 4 and, then 5 . According to lemma 2 | ri | < 
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< I a± | , Vi = l,n. Now we shall follow the algorithm a 
second time, but this time for an equation in which ri 
(according to step 5) is substituted for a±. Again, 
according to lemma 2, the new |r| written | r 2 I will have 
the propriety: | r 2 I < | ri | . We will get to |r| = 1 as 

| r | >1 and |r| <4, and if |r| ...1, following the 
algorithm once again we get |r| < |ri|, a.s.o. Hence, the 
algorhythmus has a finite number of repetitions. 

Theorem of Correctness . The former algorithm 
calculates correctly the general integer solution of the 
linear equation (2) . 

Proof 

According to lemma 5 the algorithm is finite. From 
lemma 1 it follows that the set M has a minimum, hence 
step 2 of the algorithm has meaning. When |r| = 1 it was 
shown in lemma 3 that step 3 of the algorithm calculates 
correctly the general integer solution of the respective 
equation (the equation that appears at step 3) . In lemma 
4 it is shown that if |r| ...1, by the substitutions steps 
4 and 5 introduce in the initial equation the general 
integer solution remains unchanged. That is, we pass from 
the initial equation to a linear system having the same 
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general solution as the initial equation. The variable h 
is a counter of the newly introduced variables which are 
used to successively decompose the system in systems of 
two linear equations. The variable p is a counter of the 
substitutions of variables (the relations, at a given 
moment, between certain variables) . 

When the initial equation was decomposed to |r| =1, 
we have to follow the reverse way: i.e., to compose its 

general integer solution. This reverse way is directed by 
the substeps 3 (A) , 3 (B) and 3 (C) . The substep 3 (D) has 
only an aesthetic role, i.e., to have the general solution 

under the form: x± = fi(ki, ..., k n -i) , i = l,n, f± being 

linear functions with integer number coefficients. This 
"if possible" shows that substitutions are not always 
possible. But when they are we have to make all the 
possible substitutions. 

Note 1 . The former algorithm is written under a form 
easy to introduce in the computer. 

Note 2 . The former algorithm is more "rapid" than 
that of W. Sierpinski's 1, i.e., the general integer 
solution is reached after a smaller number of iterations 
(or, at least, the same) for any linear equation (2) . In 
the first place, both aim at obtaining the coefficient + 1 
for at least one unknown variable. While Sierpinski 
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started only by chance by decomposing the greatest 
coefficient in the module (writing it under the form of a 
sum between a multiple of the following smaller 
coefficient (in the module) and the rest) , in our 
algorithm this decomposition is not accidental but always 
seeks the smallest |r| and also chooses the coefficients 
ai and aj for which this minimum is achieved. That is, we 
test from the beginning the shortest way to the general 
integer solution. Sierpinski does not attempt at finding 
the shortest way; he knows that his way will take him to 
the general integer solution of the equation and is not 
interested in how long it will be. However, when an 
algorithm is introduced in a computer, it is preferable 
that the computer time should be as short as possible. 

Example 1 

Let us solve in Z J the equation: 17x - 7y + lOz = 

= -12. We apply the former algorithm. 



1 . 


h = 


l, p 


= 1 


2 . 


r = 


3, i 


=3, j = 2 


3. 


1 3 | 


... 1, 


go on to step 4 . 


4 . 


(1) 








y = 


tl = 


10-3 

! z = ti + 

-7 
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5. Consider y: = ti h: = 2 

a3 : = 3 p : = 2 

the other coefficients and variables remaining 
unchanged, go back to step 2 . 



2. r = -1, i = 1, j = 3 

3. I -1 I = 1 

x = -l(-3t 2 - ( — 7 1 1 ) -12) = 3t 2 - 7t i - 12 

17- (-1 ) -1-17 

z = - 1 ( 1 7 1 2 + ( — 7 1 1 ) ! + (-12)) = 

3 3 

= -17t 2 + 42ti - 72 

(A) We substitute the values of x and z thus 
determined in the only relation (p) we 
have : 

(1) y = ti + z = -17t 2 + 43ti - 72 

(B) The substitution is not possible. 

(C) The substitution is not possible. 

(D) The general integer solution of the 
equation is: 

x = 3ki - 7k 2 + 12 
y = -17ki + 43k 2 - 72 
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z = -17ki + 42k 2 - 72; ki, k 2 0 Z 
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INTEGER NUMBER SOLUTIONS OF LINEAR SYSTEMS 
Definitions and Properties of the 
Integer Solution of a Linear System 



n 

Let 6 aijXj = bi, i = l,m (1) 

j=l 



A linear system with all its coefficients' integer 

numbers (the case with rational coefficients is reduced to 
the same) . 

o 

Definition 1. Xj = Xj, j = l,n is a particular 

o n o 

integer solution of (1) if Xj 0 Z, j = l,n and 6 aijXj = 

j = l 



= bi, i = l,m. Let the functions be fj : Z h 6 Z, j = l,n, 
where h 0 N* . 



Definition 2 . Xj = fj (ki, ...k h ), j = l,n, is the 
general integer solution for (1) if: 

n 

(a) 6 aijfj (ki, k h ) = bi, i = l,m 

j=l 

of (ki, . . . , k h ) 0 Z; 



(b) 



k h ) 0 



o 

For any Xj = Xj, j = l,n, particular 

o o 

integer solution of (1), there is (ki, ..., 



o 

k h ) = 



= 1 , n . 



o 

0 Z so that fj (ki. 



• • • r 



x j , Vj 
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(In other words, the general solution is 
the solution that comprises all the other 
solutions . ) 

* 



Propriety 1 

A general solution of a linear system of m equations 
with n unknowns, r (A) =m<nisn-m times 
undetermined . 

Proof 

We assume by reductio ad absurdum that it is of order 
r, 1 < r < n - m (the case r = 0, i.e., the solution is 
particular, is trivial) . It follows that the general 
solution is of the form: 

Xi = UnPl + ... + UlrPr + Vi 

(Si) 

X n = UnlPl + . . . + UnrPr + V n , U ih , V± 0 Z 

p h = parameters 0 Z 

We prove that there are n - m times undetermined 
solutions. The homogenous linear system (1), solved in r 



admits the solution: 
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Dm+l 

Xi = X m+ 1 + 

D 



D n 

+ 

D 



Xn 



• Dm+1 

Xm X m+ i + 

D 



+ 



m 




X n 

D 



o 

Let xi = xi, i = l,n be a particular solution of the 
linear system (1) . 

Considering 



x m+ i = D ! k m+ i, we get a solution 

• i i 

• xi D m+ ik m+ i + . . . + D n k n “t - 

Xn D ! k n 

• m m 

Xm Dm+lkm+1 i . . . 1 D n k n 

0 

x m+ i = D ! k m+ i + x m+ i 



• 0 

Xn D ! kn 1 X n , kj 

= parameters 0 Z 



which depends on the n - m independent parameters, for 



o 

Xi 



o 

x m 



the 



system (1) . Let the solution be n 



m times undetermined: 
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Xi Cuki T ... "i Cin-mkn-m 4* di 

(S 2 ) 

X n Cnlki + ... 4" Cnn-mkn-m 4 dn 

Cij, di 0 Z, kj = parameters 0 Z 

(There are such solutions, we have proven it before.) 
Let the system be: 



anxi + ... + ai n x n — bi 



a±j, b± 0 Z 

a m iXi + . . . + a mn x n = b m xi = unknowns 0 Z 



I. The case bi = 0, i = l,m results in a homogenous 
linear system: 

anXi + . . . + ain =0, i = l,m 

(S 2 ) i an (Cuki 4 ... "i Cin-mkn-m 4 di ) 4 ... ain (Cnlki 1 

"t ... Cnn-mkn-m 



+ d n ) 



0 
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0 — (aiiCn + ... + ainCnl)ki + ... + (ailCin-m + 

+ ... + 9-inCnn-m ) ! kn-m + (aildi + . . . + 

+ a in dn) , V kj 0 Z 

For ki = ... = k n -m = 0 Y andi + . . . + ai n d n = 0 

For ki = ... = k h -i = k h +i = . . . = k n -m = 0 and k h = 

= 1 Y (anCih + ... + ainCnh) + (audi + ... + 

Y a in d d n) 0 Y anCih y ... y ainCnh 0 

V i = l,m, V h = l,n-m 



Vect. V h = 
system . 



Cih | 



|, h = l,n-m, are the particular 

• I 

Cnh I solutions of the 



V h , h = l,n-m are also linearly independent because the 
solution is n - m times undetermined. {Vi, ..., V n - m } + d 
is a linear variety that includes the solutions of the 
system obtained from (S 2 ) . Similarly, for (Si) we deduce 



that 



s 



are particular solutions of 
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the given system and are linearly independent because 
(Si) = r - times undetermined solution and 

Vi 

V = . = a solution of the given system. 

Vn 

The case (a ) . Ui, U r , ce = linearly dependent, 

it follows that {Ui, U r } is a free submodule of order 

r < n - m of solutions of the given system, then, it 
follows that there are solutions that belong to {Vi, ..., 
Vn-ml + d and which do not belong to {Ui, ..., U r }, a fact 
which contradicts the assumption that (Si) is the general 
solution . 

The case (b ) . Ui, . .., U r , V = linearly independent. 

{Ui, ..., U r } + V is a linear variety that comprises the 
solutions of the given system, which were obtained from 
(Si) . It follows that the solution belongs to {Vi, . . . , 
Vn-m} + d and does not belong to {Ui, ..., U r } + V, a fact 
which is in contradict ion with the assumption that (Si) is 
the general solution. 



II . 



When there is an i 0 l,m, bi ... 0 Y nonhomogeneous 
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linear system 

anxi + . . . + ainX n = bi, i = l,m 
(S 2 ) Y an (Cuki + ... + Cin-mkn-m + di ) + ... + 

"t ain (Cnlki 1 ... + Cnn-mkn-m dn bi^ 

it follows that 

Y (ailCii + ... + ainCnl) ki + ... + (auCin-m + 

"t ... "t ainCnn-m) k n -m + (andi 1 ... + ainQ-n) 

bi , 

for ki = ... = k n -m = 0 Y andi + . . . + ai n d n = bi; 

for ki = ... = kj-i = kj+i = . . . = k n -m = 0 and kj = 1 Y 

Y (anCi j + ... + 3-inCnj + (a in di + . . . + 

+ ain d n ) = bi it follows that 

ailCij 1 ... + ainCnj 0 

aiidi + . . . + aindn = bi, V i = l,m, V j = l,n-m; 



Clj 
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Vj = . , j = l,n-m are linearly independent 

c n j because the solution (S 2 ) in 

n - m times is undetermined. 

di 

?! Vj, j = l,n-m and d = . are linearly independent. 

d n 

We assume that they are not linearly independent. It 
follows that 

d S 1 Vi 1 ... + Sn-mVn-m 

SlCn "t ... + Sn-mCin-m 

SlCnl I ... I Sn-mCnn-m 



Irrespective of i = l,m: 



bi andi + ... f Sindn an (SiCn + ... 1 Sn-mCin-m) 



t ... I Sin (SlCnl I ... I Sn-mCnn-m) ( SnC f ... + 



+ ainCnl) Si 1 ... 1 (anCin-m f ... + ainCnn-m) , S n -m 




50 



Then, bi = 0, irrespective of i = l,m, contradicts the 

hypothesis (that if there is an i 0 l,m, bi ... 0) . It 
follows that Vi, V n -i„, d are linearly independent. 
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{Vi, Vn-m} + d is a linear variety that contains 

the solutions of the nonhomogeneous system, solutions 
obtained from (S 2 ) • Similarly, from (Si) it follows that 
{Gi, ..., G r } + V is a linear variety containing the 
solutions of the nonhomogeneous system, obtained from 
(Si) . 

n-m > r follows that there are solutions of the 
system that belong to {Vi, . .., V n - m } + d and which do not 
belong to {Gi, ..., G r } + V (it contradicts the fact that 
(Si) is the general solution) . Then, it results that the 
general solution depends on the n-m independent 
parameters . 

Theorem 1 . The general solution of a nonhomogeneous 
linear system is equal to the general solution of an 
associated linear system plus a particular solution of the 
nonhomogeneous system. 

Proof 

Let the homogeneous linear solution: 

anXi + . . . + ainXn = 0 

, (AX = 0) 



3-mlXi + ... T UmnX n 



0 
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with the general solutions: 

Xi Cuki + ... + Cin-mkn-m di 

X n Cnlki + ... + Cnn-mkn-m 1 dn 

0 

and xi = xi 

• o 

X n = X n 

a particular solution of the nonhomogeneous linear system 
AX = b; 

o 

?! Xi = Cuki + . . . + Cin-mkn-m + d + Xi 

• 0 

Xn Cnlki + ... + Cnn-mkn-m 1" dn f X n 

is a solution of the nonhomogeneous linear system. 

We have written 

an ... ain xi bi o 

A = . , X = . , b = . , 0 = . 

a m i ... a^n Xn b m o 
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(vector of dimension m) , 



o 



Xi 



ki 



Oil • • • Cin-m 



di 



K = 



, C 



, d 



, x 



0 



k 



n-m 



X n 



Cnl 



Cnn-m 



d 



n 



AX = A (Ck + d + x°) = A (Ck + d) + AX° = b + 0 = b 



o 

We will prove that irrespective of xi = yi 



• o 

Xn = y n 

there is a particular solution of the nonhomogeneous 
system 



o 

ki = ki 0 Z 



k 



n-m 



0 

kn-m 0 Z 



with the propriety: 



xi = cnki + . 



o oo 

. + Cinkn-m + di + Xi = yi 



n Cnlki + 



0 0 0 
• Cnn-mkn-m di + X n yn 



X; 
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o 

We write y° = yi 



• o 

Yn 



o 

We demonstrate that those kj 0 Z, j = l,n-m are those 
for 

o 

which A(CX° + d) = 0 (there are such Xj 0 Z because 



xi = 0 



Xn = 0 



is a particular solution of the homogenous linear system 

and X = CK + d is a general solution of the 
nonhomogeneous 

linear system) A(CK° + d + X° - Y°) = A(CK° + d) + AX° - 

- AY 0 =0+b-b=0. 

Propriety 2 . The general solution of a homogenous 
linear system can be written under the form: 



Xi Cuki 1 ... f Cin-mkn-m 

(SG) 

( 2 ) 

X n Cnlki + ... 1 Cnn-mkn-m r 
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kj = a parameter belonging to Z (with di = ch = ... = 

= d n = 0) 

Proof 

(SG) = general solution. It results that (SG) is 
(n-m) times undetermined. 

Let (SG) of the form be 

Xi CuPi f ... 1 Cin-mPn-m t di 

X n CnlPl ... “t“ Cnn-mPn-m dn 

with not all d± = 0 ; we demonstrate that it can be written 
under the form (2); the system admits the trivial solution 

xi = 0 0 Z 

• f 

x n — 0 0 Z 



it results that there are pj 0 Z, j = l,n-m 

o o 

Xi CuPi + ... + Cin-mPn-m di 0 
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(4) 

• o o 

X n CnlPl ... "I" Cnn-mPn-m d n 0 



0 

Substituting pj = kj + p jf j = l,n-m in (3) 



kj 0 Z 

o 

Pj 0 Z A Y pj 0 Z 



Pj 0 Z 

o o 

Pj 0 Z A Y kj = pj - pj 0 Z 



which means that they do not make any restrictions. 
It results that 



o 

Xi = Cuki + ... + Cin-mkn-m + (CnPi + 



0 

Cin-mPn-m di ) 



0 

X n C n iki + ... + Cnn-mkn-m ( C n lPl 



+ 



0 

Cnn-mPn-m 



+ d n ) 



0 0 

But ChlPl + ... + Chn-mPn-m + d h = 0, h = 1,11 (from (4)) . 

Then the general solution is of the form: 



Xi 



Cuki + ... + Cin-mkn-m 
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Xn Cnl^i + ... + Cnn-mkn-m 



kj = parameters 0 Z, j = l,n-m; it results that di = d 2 = 
= ... = d n = 0 

Theorem 2 . Let the homogenous linear system be: 

anXi + . . . + ainX n = 0 

a m ixi + ... + a mn x n = 0, r (A) = m; 

(a h i, . . . , a hn ) = 1, h = l,m and the general solution 

Xi Cuki + ... 1 Cin-mkn-m 

Xn Cnlki + ... + Cnn-mkn-m 

then ( Uhl r . . . r Uhi-lf 3-hi + lf • • • r Shn) I (Cii^ . . . t Cin-m) 

irrespective of h = l,m and i = l,n. 

Proof 

Let some arbitrary be h 0 l,m and some arbitrary i 0 

1 t n/ a h lXi + ... + ahi-lXi-i + 3-hi + lXi + i + ... + ShnXn ahiXi . 
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Because (ahi, . . . , ahi-i, ahi+i, . . . , ahn) I ahi it results that 
d = (ahi, . .., ahi-i, ahi+i, •••, ahn) |xi, irrespective of 
the value of xi in the vector of particular solutions; for 
k 2 = k 3 = ... = k n -m = 0 and ki = 1 we get the particular 
solution : 

Xi = Cii 

xi = Cii Y d | On and so on 

Xn Cnl 

for xi = k 2 = ... = kn-m-i = 0 and k n - m = 1 the following 
particular solution results: 

Xi Cin-m 
Xi Cin-m 
Xn Cnn-m 



it results that d|Ci n - m ; hence, d|Cij, j = l,n-m Y 
Y d | ( Cii , . . . , Cin-m) • 
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Theorem 3 . 

xi = Cuki + . . . + Cin-mkn-m Cij 0 Z being given 

If 

x n = Cniki + . . . + Cnn-mkn-m, kj = parameters 0 Z 
is the general solution of the homogenous linear system 
anXi + . . . + ainXn = 0 m < n; 

a m ixi + . . . + a mn x n = 0 r (A) = m; 



then (Cij, c n j) = 1, V j = l,n-m. 



Proof 



We assume, by reductio ad apsurdum, that there is 



jo 0 l,n-m: (ch , c n i ) = d ... 1, we consider the 

o o 

maximal co-divisor >0; we reduce the case when the 

maximal co-divisor is -d to the case when it is equal to d 
(nonrestrict ive hypothesis) ; then the general solution can 

be written under the form: 



Cuki + ... i Cij dkj i ... + Cin-mkn-m 

0 0 



Xi 
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( 5 ) 



X n = Cnlki + 



+ Cni dki + 
0 0 



"t Cnn-mkn-m 



where d = (cii , . . . , c n i ) , ch = d ! cn and (cii 

o o o o c 

• • • r Cni ) 1 • 

0 

We prove that 



Xi = Cij 

o 



Xn Cn j 

0 



is a particular solution of the homogenous linear system. 
We write 



Cll ... Clj d ... Cin-rr 



c = . 



k = 



Cni ... Cnj d ... Cnn-rr 



kn 



k n -rr 



x = c ! k the general solution; 



an ... ai 
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We know AX = 0 Y A(CK) = 0, A = . ; 

3-nl ... 3-mn 

We assume that the principal variants are xi, . .., x m (if 
not we have to renumber) . It follows that x m+ i, ... x n is 
the secondary variant. 

For ki = ... = k-j -i = k-i +i = ... = kn-m = 0 and kj = 1 

o o o 

we get a particular solution of the system 

xi = Ci-j d Ci-j d 

o 

• • 0 

Y 0 = A . = d ! A 

Xn Cnj d Cnj d 

0 

0 

? lj ° 

Y A . = 0 Y 

Cnj 

0 

is the particular solution of the system. 

We demonstrate that this particular solution cannot 
be obtained by 

Xi Cuki f ... Cij dkj f ... + Cin-mkn-m Cij 

0 0 



Xi = Cij 

• 0 



Xn Cnj 

0 




( 6 ) 



o 
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Xn Cnl^i + ... + Cnj dkj + ... + Cnn-m^n-m Cnj 

0 0 0 



Xm+1 Cm+lki + 



1 Cm+idki + ... + Cm+1 . n-mkn-m 

0 



Cm+l^j 



(7) 



X n Cnlki + ... + Cni dk-i ... + Cnn-mkn-m Cnj 



■'HI 

0 0 



Cm+1, 1 ... Cm+lj 



Cm+ln-rr 



Chi • • Cnj ... Cn,n-m 1 

Y kj = = - 6 Z (because 

d...l ) . 

0 . d 

Cm+1,1 ... Cm+l,j d ... Cm+l,n-m 
0 



Cnl 



Cn 



d 



Cn 



0 

It is important to point out the fact that those kj = 



j = l,n-m that satisfy system (7) also satisfy system 

( 6 ), 

because, otherwise (6) would not satisfy the definition 

of the solution of a linear system of equations (i.e., 

considering system (7) the hypothesis was not 

restrictive) . From Xj 0 Z it follows that (6) is not 
the 

o 

general solution of the homogenous linear system 
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contrary to the hypothesis) ; then (cij . . . c n j) = 1, 

irrespective of j = l,n-m. 

Propriety 3 . Let the linear system be 

anxi + ... ai n x n = bi 

S-mlXi I ... "t amnXn bm 

aij, bi 0 Z, r (A) = m < n, Xj = unknowns 0 Z 
Solved in R, we get 

Xi fl (X m+ i, . . . t X n ) 

; xi, . . . , x m are the 

Xm fm (X m+ i, . . . t X n ) 

main variants where f± are linear functions of the form: 



Cm+lXm+d ... +CnX n “t“Si i 

f± = where c m+ j, di, ei 0 Z; i = l,m; 

di 



j = l,n-m. 



Si 

-- 0 Z irrespective of i = l,m then the linear system 



If 
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di 

admits integer solution. 

Proof 

For 1 < i < m, x± 0 Z, then fj 0 Z. Let: 

Xm+1 Um+lkm+1 



X n 

Solution 



Xi 






= U n k n 



1 1 ei 

v m+ ik m+ i + ... f v n k n F 

di 

m m Sm 

V m+ ik m+ i 4" ... + v n k n "t 

dm 



where u m +i is the maximal co-divisor of the denominators 
of 

i 

Cm+j 

the fractions , i = l,m, j = l,n-m calculated after 

di 

their simplification when they were irreducible. 



i Cm+jUm+j 

v m+ j = 0 Z; this is a solution n-m times 

di 



undetermined (depends on n-m independent parameters: 



k m+n , . .., k n ) but is not a general solution. 
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Propriety 4. Under the conditions of propriety 3, 
if 



io 

there is an i 0 0 l,m: = u m+ ix m+ i + 

o 



ei 

io 

+ u n x n + 

O: 

0 
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e± 

io 

with u m+ j 0 Z, j = l,n-m and 6 Z then the system 

does 

di 

o 

not admit integer solution. 



Proof 

ce x m+ i, . .., x n in Z it results in xi 6 Z. 

o 

Theorem 4 . Let the linear system be 



anxi + ... + ai n x n — bi 
9-mlXi + ... + amnXn bm 

aij, b± 0 Z, Xj = unknowns 0 Z, r (A) = m < n. If there 
are 

indices 1 < ii < ... < i m < n, i h 0 {1, 2, ..., n}, 

h = l,m with the propriety: 



aii 

i 



au 

m 



and 



3-mi 

1 



3-mi 

m 
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a 

x 



ii 



bi 



aii 

2 



bm 



ami 



2 



an 

m 

is divided by 

ami 

m 



. aii 

a 1 

x 

lm 



an bi 

m-1 

is divided by 



ami 



1 



ami 

m-1 



bm 



then the system admits integer number solutions. 



Proof 

We use propriety 3 



d± 



l,m; e± = a x , h = l,m 

h 

lh 



Note 1 . Conversely, it is not true. 

Consequence 1 . Any homogenous linear system admits 
integer number solutions (beside the trivial one); r (A) 

= m < n . 



Proof 
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a x = 0 ! irrespective of h = l,m. 
ih 

Consequence 2 . If - = + 1, it follows that the 
linear system admits integer number solutions. 

Proof 



! (+1), irrespective of h = l,m; 

ih 



s x 0 Z. 
ih 
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FIVE INTEGER NUMBER ALGORITHMS TO 
SOLVE LINEAR SYSTEMS 

This chapter further extends the results obtained in 
4 and 5 (from linear equations to linear systems) . Each 
algorithm is strictly demonstrated and then an example is 
given . 

Five integer number algorithms to solve linear 
systems are further given. 

Algorithm 1 (method of substitution) 

(Although simple, this algorithm requires complex 
calculus but is, nevertheless, advantageous in introducing 
it in the computer) . 

Some integer number equations are initially solved 
(which is usually simpler) by means of one of the 
algorithms 4 or 5 . (If there is an equation of the system 
which does not admit integer systems, then the integer 
system does not admit integer systems. Stop.) the 
general integer solution of the equation will depend on n- 
1 integer number parameters (see 5) : 

(i) (i) (i) 

(Pi) xi = fi (ki, ..., k n -i ), i = l,n where all 

i i 
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(i) 

the functions fi are linear with integer number 
coefficients . 

This general integer number system (pi) is introduced 
in the other m-1 equations of the system. We get a new 
system of m-1 equations with n-1 unknown variables: 

(i) 

ki , ii = l,n-l, which is also to be solved in integer 
numbers. The procedure is similar. Solving a new 
equation, we obtain its general integer solution: 

( 1 ) ( 2 ) ( 2 ) ( 2 ) 

(p 2 ) ki = fi (ki , ..., k n -2 ), i 2 = 1 / n-1 

2 2 



( 2 ) 

where all the functions fi are linear, with integer 

2 

number coefficients. (If, along this algorithm we come 
across an equation which does not admit integer solutions, 
then, the initial system does not admit integer solution. 
Stop . ) 

In the case that all the solved equations admitted 
integer systems at step (j), 1 < j < r, (r being of the 

same rank as the matrix associated to the system) then: 

(j-D (j) (j) (j) 
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(p j ) ki = f± (ki , k n -j ) , ij = 1 , n- j + 1 , 

j j 

(j) 

fi are linear functions with integer number 

j 

coefficients . 

Finally, after r steps, and if all the solved 
equations admitted integer solutions, we get to the 
integer solution of one equation with n-r+1 unknown 
variables . 

The system will accept integer solutions if and only 
in this last equation will have integer solutions. If it 
does, let the general integer solution of it be: 



(r-l) (r) (r) (r) 

(p r ) ki = fi (ki , ..., kn-i) , ir = 1 , n-r+1. 



(r) 

where all fi are linear functions with integer number 

r 

coefficients . 

Now the reverse way follows. 



( r-l ) 

We introduce the values of ki , i r = 1, n-r+1 



( r-2 ) 

at step (p r ) in the values of ki , i r _i = l,n-r+2 from 

(r-l) 

step (p r -i) . 



It follows: 
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( r-2 ) ( r-1 ) (r) (r) (r) (r) (r) (r) 

k-i f i ( f 1 (ki ^ f kn-r^ • > • r fn-r+1 (kl f / kn-r) ) 

r-1 r-1 



(r-1) (r) (r) 

CJi (ki r • • • r kn-r ) r 1 r - ' 1 r H 1 1 

r-1 



(r-1) 

from which it follows that gi are linear functions with 



integer number coefficients. 

Then follow those from (p r - 2 ) in (p r - e ) , and so on, 
until we introduce the values obtained at step (P 2 ) in 
those from the step (pi) . It will follow: 



(1) (r) (r) 

xi = gi (ki , . . . , k n -r notation gi (ki, . . . , k n - r ) , 



i = 1 , n 

with all gi , most obviously, linear functions with 
integer number coefficients (the notation was made for 
simplicity and aesthetical aspects) . This is, thus, the 
general integer solution, of the initial system. 

The correctness of algorithm 1 . The algorithm is 
finite because it has r steps on the first way and r-1 
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steps on the reverse. (r < + 4) . Obviously, if one 
equation of one system does not accept (integer number) 
solutions then the system does not have either. 

Writing S for the initial system and Sj the system 
resulted from step (p j ) , 1 < j < r-2, it follows that 
passing from (pj) to (pj+i) we pass from a system Sj to a 
system Sj+i equivalent from the viewpoint of the integer 



(j-l) o 

number solution, i.e., ky = t± , ij = l,n-j+l which is 

j j 

a particular integer solution of the system Sj if and only 

(j) o 

if ki = hi , ij+i = l,n-j is a particular integer 

j+i j+i 

0 (j+l) 0 0 

solution of the system Sj+i where hi = f± (ti,...,t n - 



ij+i = l,n-j. Hence, their general integer solutions 
are 

also equivalent (considering these substitutions) . So 

that, in the end, the solving of the initial system S is 

equivalent with the solving of the equation (of the 
system 

consisting of one equation) S r -i with integer number 

coefficients. It follows that the system S admits 
integer 

number solution if and only if all the systems Sj admit 
integer number solution, 1 < j < r-1 . 

Example 1 . By means of algorhythmus 1, let us 
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calculate the integer number solution of the system: 

5x - 7y - 2 z + 6 w = 6 
(S) 

9 -4x + 6 y - 3z + llw = 0 

Solution: We solve the first integer number equation. 

We 

obtain the general integer solution (see [4] or [5]): 

x = ti + 2 t 2 
y = ti 

(Pi) 

z = -ti + 5t 2 + 3 t 3 - 3 
w = t 3 , 

where ti, t 2 , t 3 0 Z. 

Substituting in the second, we get the system: 

(Si) 5t i — 23t 2 + 2t 3 + 9 = 0 

Solving this integer equation we obtain its general 
integer solution: 

ti = ki 

(p 2 ) t 2 = ki + 2 k 2 + 1 

t 3 = 



9ki + 23k 2 + 7 
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where ki, k 2 0 Z . 

The reverse way . Substituting (p 2 ) in (pi) we 
obtain : 

x = 3ki + 4k 2 + 2 
y = ki 

z = 31ki + 7 9k 2 + 23 
w = 9ki + 23k 2 + 7 , 

where ki, k 2 0 Z which is the general integer solution 
of 

the initial system (S) . Stop. 

Algorithm 2 
Input 

A linear system (1) without all aij = 0. 



Output 

We decide on the possibility of an integer solution 

of this system. If it is possible, we obtain its 
general 

integer solution. 

Method 

1. t = 1 , h = 1 , p = 1 
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2 . 

unknown 



the 



(A) Divide each equation by the maximal co- 
divisor of the coefficients of the 

variables. If you do not get an integer 
quotient for at least one equation, then 

system does not admit integer solutions. 



Stop . 

(B) If there is an inequality in the system, 
then the system does not admit integer 
solutions. Stop. 

(C) If the repetition of more equations 

occurs , 

keep one and if an equation is an 

identity, 

remove it from the system. 

3. If there is (i 0 , j 0 ) so that | a± i | =1, then 

oo 

obtain the value of the variable from the 

o 

equation i 0 ; relation (T t ) . Substitute this 
relation (where possible) in the other 

equations 

of the system and in the relations (T t -i) , (H h ) 
and (P p ) for all i, h and p. Consider t: = t 

+ 1 ' 



there 



remove equation (i 0 ) from the system. If 



is not such a pair, go on to step 5. 
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4. Does the system (left) have at least one 
unknown 

variable? If it does, consider the new date 

and 

go on to step 2. If it does not, write the 
general integer solution of the system 
substituting ki; k 2 , . . . for all the variables 

from the right term of each expression which 
gives the value of the unknowns of the initial 
system. Stop. 

5. Calculate a = min {|r|, an / r (mod an ), 

1 2 

h Jl, 3 2 

0 < | r | < | a±i | }, and determine the indices 

2 

i, ji, j 2 as well as the r for which this 

minimum 

can be calculated. (If there are more 

variants , 

choose one, arbitrarily.) 

a±j -r 

i 

6. Write: Xj = th xij , relation (H h ) . 

2 = . . 1 
dig 

2 

Substitute this relation (where possible) in 

all 

the equations of the system and in the 

relations 

(T t ) , (H h ) and (P p ) for all t, h and p. 

7. (A) If a ... 1, consider Xj : = t h , h: = h + 1 
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and 

2 

go on to step 2 . 

(B) If a = 1, then obtain the value of Xj 

from 

i 

from the equation (i) ; relation (P p ) . 
Substitute this relation (where possible) in 

the 

other equations of the system and in the 
relations (T t ) , (H h ) and (P p -i) for all t, h 

and p . 

Remove the equation (i) from the system. 
Consider h : = h + 1 , p : = p + 1 and go back to 
step 4 . 

The correctness of algorithm 2 . Let the system (1) 

be . 

Lemma 1 . We consider the algorithm at step 5. 

Also, 

let M = { | r | , a ±j / r (mod a±j ) , 0 < |r| < | aij |, i. 



j 2 = 1 , 2, 3, . . . } . Then M ... . 

Proof 

Obviously, M is finite and M d N* . Then, M has a 
minimum if and only if M ... . We suppose, conversely, 
that M = . Then a±j / 0 (mod a±j ), Vi, j x , j 2 . It 

follows conversely as well that a±i / 0 (mod a±j ) , V i, 

2 1 

ji, j 2 . That is | aij | = I a ±j | , V i, j x , j 2 . We 
consider 

1 2 

a io arbitrary but fixed. It is clear that (a± i, . . ., 

o 
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. . a± n ) - a± j ... 0, V j (because the algorithm has 
passed 

o o 

through the substeps 2 (B) and 2 (C) ) . But, as it has 
also 

passed through step 3, it follows that | ai j | ... 1, V j, 
but 

o 

as it previously passed through step 2 (A) , it would 
result 

that | ai j | = 1, V j. This contradict ion shows that the 
assumption is false. 

Lemma 2 . Let a± / r (mod a±i ) . Substitute Xj = 

0 1 2 2 

a± j -r 

0 1 

= th Xj in the system (A) obtaining the system 



o 

(B) . Then, Xj = Xj, j = l,n is the particular integer 

o 

solution of (A) if and only if Xj = Xj, j ... j 2 and th = 



o 




ai j -r 

0 1 



o 2 



+ 



is the particular integer solution of 
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Lemma 3 . Let ai ... 1 and a 2 be obtained at step 5. 



Then 0 < a 2 < ai . 



Proof 

It is sufficient to show that ai < | aij I , V i, j 
because 

in order to get a 2 step 6 is obligatory, when the 
coefficients of the new system are calculated, a 2 being 
equal to a coefficient from the new system (equality of 
modules) , the coefficient on (i 0 ji) . 

Let a± j with the propriety | a± j | < ai . Hence, ai 

> 

— 0 0 0 0 

> | a± j | = min { I ay j I } . Let ai j with | a± j | > | a±j | 
there 

0 0 0 s 0 s m 

is such an element because | a± -i | is the minimum of the 

0 m 

coefficients in the module and not all | a± j | , j = l,n 
are 

o 

equal (conversely, it would result that (ai j, . . ., ai n ) 

o o 

- a± j, V j 0 l,n; the algorithm passed through substep 
2 (A) 

o 

has simplified each equation by the maximal co-divisor 
of 

its coefficients; hence, it would follow that | ai i I = 1, 

o 

V j = l,n which, again, cannot be real because 

the algorithm has also passed through step 3) . Of the 




81 



coefficients a± j we choose one with the propriety ai j 

0 m Os 

0 



... ;a± j there is such an element (contrary, it would 
result 

0 m 

(a± i, ... ain) - I a± j | , but the algorithm has also 
passed 

0 0 0 m 

through step 2 (A) and it would mean that | a± j | =1 
which 

0 m 

contradicts step 3 through which the algorithm has also 
passed) . 

Considering q 0 = a± j / a± j 0 Z and r = ai j 

Os 0 m Os 

" 0 _ 0 

- q 0 ai i 0 Z, we have ai -j / r 0 (mod a± -j ) and 0 < 

0 m Os 0 m 

0 

< | r 0 1 < I ai j | < | ai j I < ai . We have, thus, obtained 
an 

0 m 0 0 

r 0 with | r 0 I < ai, which is in contradiction with the 
very 

definition of ai . Thus, ai < | aij | , V i, j . 

Lemma 4 . Algorithm 2 is finite. 



Proof 

The functioning of the algorithm is meant to 
transform a linear system of m equations and n unknowns 
into one of mix ni, with mi < m, ni < n and, thus, 
successively into a final linear equation with n - r + 1 
unknowns (where r is the rank of the associated matrix) . 
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This equation is solved by means of the same algorithm 
(which works as [5]) . The general integer solution of the 
system will depend on the n - 1 integer number independent 
parameters (see [6] --similar proprieties can be 
established also for the general integer solution of the 
linear system) . The reduction of equations occurs at 
steps 2, 3 and substep 7(B) . Steps 2 and 3 are obvious 
and, hence, trivial; they can reduce the equations of the 
system (or even put an end to it) but only under 
particular conditions. The most important case finds its 
solution at step 7 (B) , which always reduces one equation 
of the system. As the number of equations is finite, we 
come to solve a single integer number equation. We also 
have to show that the transfer from one system itux n± to 
another m i+ ixni+i is made in a finite interval of time: by 

steps 5 and 6 permanent substitution of variables are made 
until we get to a = 1 (we get to a = 1 because, according 
to lemma 3, all a - s are positive integer numbers and 
form a strictly decreasing row) . 

Theorem of correctness . Algorithm 2 correctly 
calculates the general integer solution of the linear 
system . 



Proof 




83 



Algorithm 2 is finite according to lemma 4. Steps 2 
and 3 are obvious (see also [4], [5]). Their part is to 

simplify calculus as much as possible. Step 4 tests the 
finality of the algorithm; the substitution with the 
parameters kg, k 2 , . . . has systemization and aesthetic 

reasons. The variables t, h, p are counter variables 
(started at step 1) and they are meant to count the 
relations of the type T, H, P (numbering required by the 
substitutions at steps 3, 6 and substep 7 (B) ; h also 
counts the new (auxiliary) variables introduced in the 
moment of decomposition of the system. The substitution 
from step 6 does not affect the general integer solution 
of the system (it follows from lemma 2) . Lemma 1 shows 

that at step 5 there is always a, because ... M d N* . 

The algorithm performs the transformation of a system 
rriixni into another, ai+ixni+i, equivalent to it, preserving 
the general solution (taking into account, however, the 
substitutions) (see also lemma 2) . 

Example 2 . Calculate the integer solution of: 

- 12x - 7y + 9z =12 

- 5y + 8z + lOw = 0 

Oz + Ow = 0 

15x + 21z + 69w = 3 
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Solution 

We apply algorithm 2 (we purposely looked for an 
example to be passed through all the steps of this 
algorithm : 

1. t = 1 , h = 1, p = 1 

2. (A) The fourth equation becomes: 5x + 7z + 

2 3w = 1 

(B) — 

(C) Equation 3 is removed. 

3. No; go on to step 5. 

5. a = 2 and i = 1, ji = 2, j 2 = 3 and r = 2. 

6. z = ti + y, the relation (Hi) . Substituting it 

in the system: 

- 12x + 2y + 9ti = 12 

3y + 8ti + lOw = 0 

5x + 7y + 7ti + 23w = 1 

7. a ... 1; consider z: = ti, h: = 2 and go back to 

step 2 . 



2 . 

3. No. Step 5. 

5. a = 1 and i = 2, ji 



= 1 . 



= 4, j 2 = 2 and r 
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6. y = t 2 - 3w, the relation (H 2 ) . Substituting in 
the system: 



12x + 


2 1 2 


+ 


ft 

1 


6w = 


12 




3t 2 


+ 


+ 

-p 

00 


w = 


0 


5x + 


7t 2 


+ 


7ti + 


2w = 


1 



Substituting it in relation to (Hi), we get: 
z = ti + t 2 - 3w, relation (Hi) ' . 

7. w = - 3t 2 - 8ti, relation (Pi) . 

Substituting it in the system, we get: 

- 12x + 20t 2 + 57ti = 12 
5x + t 2 - 9ti = 1 

Substituting it in the other relations, we get: 



10t 2 


+ 2 5t i 


, (Hi) ’ 


1 » . 
f 


10t 2 


+ 24ti 


, (H 2 ) ' 


1 » . 
f 



h: = 3, p: = 2 and go back to step 4. 

4 . Yes 

2 . 

3. t 2 = 1 - 5x + 9ti, relation (Ti) . 

Substituting it (where possible) we get: 
{- 112x + 237ti = -8 (the new system); 

z = 10 - 5 Ox + 1 1 5t i , (Hi) ' ' ' 
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10 


- 50x + 


1 1 4t i , 


(H 2 


-3 


+ 1 5x - 


35ti, 


(Pi 



Consider t: =2; go on to step 4. 



4. Yes. Go back to step 2. (From now on algorithm 
2 works similarly with that from [5], with the 
only difference that the substitutions must also 
be made in the relations obtained up to this 
point) . 

2 . 

3. No. Go on to step 5. 

5. a = 13 (one three) and i = 1, ji = 2, j 2 = 1 and 
r = 13 . 

6. x = t 3 + 2ti, relation (H 3 ) . 

After substitution we get: 



-1 12t 3 + 13ti = -8 
z = 10 - 50t 3 + 15ti , 

y = 10 - 50t 3 + 14ti , 

w = -3 + 15t 3 - 5ti , 

t 2 = 1 - 5t 3 - ti , 



(the system) 
(Hi) IV ; 

(H 2 ) ' ' ' 

(Pi) ' ' 

(Ti) ' 



7. x: = t 3 , h: = 4 and go on to step 2. 

2 . 

3. No, go on to step 5. 
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5 . 


a = 


= 5 and 


i = 1 / 


ji = 1 / 


j 2 = 


2 and r = 


5 


6 . 


ti 


= 1 4 + 


9t 3 , relation 


(H 4 ) . 








Substituting it. 


we get 


: 5t 3 


+ 1 3 1 4 = ■ 


-8 (the 




system) . 














z 


= 10 + 


8 5 1 3 + 


15t 4 / 


(Hi) 


V 

r 






y 


= 10 + 


7 6 t 3 + 


1 4t 4 , 


(H 2 ) 


IV 

r 






X 


= 


1 9t 3 + 


2t 4 , 


(h 3 ) 


T 

r 






w 


= -3 - 


30t 3 - 


5t 4 , 


(Pi) 


! 1 1 . 
r 






t 2 


= 1 - 


1 4t 3 - 


1 4 , 


(Ti) 


T T 

r 




7 . 


ti : 


t 4 , 


h: = 5 


and go 


back to step 2 . 




2 . 


— 














3. 


No; 


step 


5 . 










5 . 


a = 


= 2 and 


i = 1 / 


ji = 2, 


32 = 


1 and r / 


-2 . 


6 . t 3 

substitution. 


= t 5 - 


3 t 4 , relation 


(H 5 ) • 


After 






we 


get : 














5t 5 


- 2 t 4 


= -8 






(the system) 




z 


= 10 


+ 8 5 1 5 


- 2 4 Ot 4 


r 


(Hi) VI ; 






y 


= 10 


+ 7 6 t 5 


- 2 1 4t 4 


r 


(H 2 ) v ; 






X 


= 


1 9t 5 


55t4 


r 


(H 3 ) iv ; 






w 


= -3 


- 30t 5 


+ 8 5 1 4 


r 


(Pi) IV ; 






t 2 


= -1 


- 1 4t 5 


+ 4 It 4 


r 


(Ti) " ' ; 






tl 


= 


9t 5 


+ 2 6 t 4 


r 


(H 4 ) ' ; 
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7 . 


t 3 : 


= t 6 , h: = 6 and go back 


to step 


2 . 


2 . 


— 








3. 


No; 


step 5 . 






5 . 


a = 


1 and i = 1 , ji = 1 , j 2 . 


r = 1 . 




6 . 


t 4 = 


= t 6 + 2t 5 , relation (H 6 ) . 


After 




substitution. 










we 


get : 








t 5 ' 


00 

1 

II 

-P 

CM 

1 


(the system) 




z 


= 10 - 395t5 - 240t6 , 


( H 2 ) VI1 


r 




y 


= 10 - 3 92t 5 - 214t 6 , 


(h 2 ) vi 


r 




X 


= - 91t 5 - 55t 6 , 


(h 3 ) ' ' 


T 

r 




w 


= — 3 + 1 4 Ot 5 + 85t6 , 


(Pl) v 


r 




t 2 


= 1 + 68 t 5 + 4 It 6 , 


(Ti) IV 


r 




tl 


= - 43t5 - 26t6 , 


(H 4 ) ' ' 


r 




t 3 


= - 5t 5 - 3t 6 , 


(H 5 ) ' 


r 


7 . 


ts = 


= 2t 6 - 8, relation (P 2 ) . 


Substituting 


in 


the 


system, we get : 0 = 0 . 







Substituting it in the other relations, it 
follows : 

z = -1030t6 + 3170 

y = - 918t 6 + 2826 
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x = - 237t6 + 728 

w = 365t6 - 1123 

1 2 = 177t6 ~ 543 

ti = 112t 6 + 344 

importance 

t3 = 1 3t 6 + 40 

1 4 = 5t6 “ 16 

Consider h: = 7, p: = 



relations of no 



3 and go back to step 4 . 



t 6 0 Z 



4. No. The general integer solution of the 
system 

is : 

x = -237ki + 728 

y = -918ki + 2826 
z = 1030ki + 3170 
w = 365ki - 1123 

where ki is an integer number parameter. 
Stop . 

Algorithm 3 
Input 

A linear system (1) . 



Output 
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We decide on the possibility of an integer solution 
of 

this system. If it is possible, we obtain its general 
integer solution. 

Method 

1. Solve the system in R n . If it does not have 

solutions in R n , it does not have solutions in 
Z n 

either. Stop. 

2. f = 1 , t = 1 , h = 1 , g = 1 

3. Write the value of each main variable x ± under 
the form: 

(Ef,i) i* Xi 0 Pij Xj "i qi + (0 rijXj + ri) /“i, 

j j 

with all q±j, qi, r ijf r±, ~i in Z so that all 

I r ±j | < \- ± \ , m ... 0, | r ± | < | | (where all Xj 
of 

the right term are integer number variables: 
either of the secondary variables of the 

system 

or other new variables introduced with the 
algorhythmus ) . For all i, we write rij / ~± . 

4. (F f ,i) i : 6 rijXj - r if j y f ,i + ri = 0 where (Y f ,i)i 

j 
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are auxiliary integer number variables. We 
remove all the equations (F f ,i) which are 
identities . 

5. Does at least one equation (F ffi ) exist? If it 

does not, write the general integer solution 
of 

the system substituting kg, k 2 , . . . for all 

the 

variables from the right term of each 

expression 

representing the value of the initial unknowns 
of the system. Stop. 

6. (A) Divide each equation (F ffi ) by the maximal 

co-divisor of the coefficients of their 
unknowns. If the quotient is not an 
integer number for at least one io then 

the 

system does not admit integer solutions. 
Stop . 

(B) Simplify--as in m--all the fractions from 
the relations (E f/i )i. 

7. Does it exist r± with the module 1? 

oo 

If it does not, go on to step 8. 

If it does, find the value of Xj from the 
equation (F ffi ) ; write (T t ) for this relation 

and 
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o 

substitute it (where it is possible) in the 
relations (E ffi ) , ( T t-1), (H h ) , (G g ) for all i. 



h and g. Remove the equation (F f ,i ) . 

Consider 

o 

f: = f + 1, t: = t + 1 and go back to step 3. 



8 . Calculate a = min { | r | , r±i / r (mod r ±i ) , 

1 2 

1 , Jl, J2 

0 < | r | < | ri-j | } and determine the indices i m , 

2 

ji, j 2 as well as the r for which this minimum 
can be obtained. (When there are more 

variants , 

choose only one) . 



a 



9. 



(A) Write Xj = z h 

2 



a± j -r 

ml 



X 



J ' 

1 



3 b 

m 2 



where z h is 



10 . 



new integer variable; relation (H h ) . 

(B) Substitute the letter (where possible) in 
the relations (E f/i ) , (F ffi ) , (T t ) , (H h _i) , 
(G g ) for all i, t, h and g. 

(C) Consider h: = h + 1. 

(A) If a ... 1, go back to step 4. 

(B) If a = 1, calculate the value of the 
variable Xj from the equation (F f ,i) ; 



i 
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relation (G g ) . Substitute it (where 
possible) in the relations (E f ,i) , (T t ) , 

(H h ), 

( G g — i ) for all i, t, h and g. Remove the 
equation (F ffi ) . Consider g: = g + 1, 
f : = f + 1 and go back to step 3 . 



The correctness of algorithm 3 

n 2 

Lemma 5 . Let i be fixed. Then ( 6 r ±j x j + ri) / s i 

j=ni 

(with all r ±j , r ir - ir ni, n 2 being integers, n 2 < n 2 , -± ... 
0 

and all Xj being integer variables) can have integer 
values 

if and only if (r in , r in , -±)/r ± . 

1 2 



Proof 

The fraction from the enunciation can have integer 

n 2 

values if and only if there is a z 0 Z so that ( 6 

rijXj + 

j=ni 

n 2 

+ ri) /~± = z ] 6 rijXj - ~i z + ri = 0 which is a linear 

j=ni 

equation. This equation admits integer solution 
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Lemma 6 . The algorithm is finite. It is true. The 
algorithm can stop at steps 1, 5 or substep 6(A) . (It 

rarely happens to stop at step 1) . An equation after 
another are gradually eliminated at step 4 and especially 
7 and 10(B) (F f ,i)— the number of equations is finite. If 

at steps 4 and 7 the elimination of equation may occur 
only in special cases, elimination of equations at the 
substep 10 (B) is always true because, through steps 8 and 
9 we get to a = 1 (see [5] ) or even lemma 4 (from the 
correctness of algorithm 2) . Hence, the algorithm is 
finite . 

Theorem of Correctness . The algorithm 3 correctly 
calculates the general integer solution of the system (1) . 

Proof 

The algorithm is finite according to lemma 6. It is 
obvious that if the system does not have solutions in R n 
it does not have in Z n either, because Z n d R n (step 1) . 

The variables f, t, h, g are counter variables and are 
meant to number the relations of the type E, F, t, H and 
G, respectively. They are used to distinguish between the 
relations and make the necessary substitutions (step 2) . 
Step 3 is obvious. All the coefficients of the unknowns 
being integers, each main variable x ± will be written: 
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Xi (0 CijXj + Ci) / i , 

j 



which can assume the form and conditions required in this 
step. Step 4 is obtained from 3 by writing each fraction 
equal to an integer variable y f ,i (this being xi 0 Z) . 

Step 5 is very close to the end. If there is no fraction 
among all (E f ,i) it means that all the main variables xi 
already have values in Z, while the secondary variables of 
the system can be arbitrary in Z, or can be obtained from 
the relations T, H or G (but these have only integer 
expressions because of their definition and because only 
integer substitutions are made) . The second assertion of 
this step is meant to systematize the parameters and 
renumber; it could be left out but aesthetic reasons 
dictate its presence. According to lemma 5 the step 6(A) 
is correct. (If a fraction depending on certain 
parameters (integer variables) cannot have values in Z, 
then the main variable which has in the value of its 
expression such a friction cannot have values in Z either; 
hence, the system does not admit integer systems) . This 
6(A) also has a simplifying role. The correctness of step 
7, trivial as it is, also results from [4], and the steps 
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8-10 from [5] or even from algorhythmus 2 (lemma 4) . 

The initial system is equivalent to the "system" from 
step 3 (in fact, (E ffi ) as well, can be considered a 
system) . So, the general integer solution is preserved 
(the changes of variables do not prejudice it (see [4], 

[5], and also lemma 2 from the correctness of algorithm 
2)) . From a system mixni we form an equivalent system 
m i+ ixni+i with m i+ i < mi and ni+i < ni. This algorithm works 
similarly to algorithm 2 . 

Example 3 . Employing algorithm 3, find an integer 
solution of the following system: 

3xi + 4x 2 + 22x 4 - 8 x 5 

6 x 1 + + 46x 4 - 12x 5 

4x 2 + 3x 3 - x 4 + 9x 5 

Solution 

1. Common solving in R 3 , it follows: 



= 25 
= 2 
= 26 



23x 4 -6x 5 -1 

Xl = 

-3 

x 4 + 2x 5 +24 

x 2 = 

4 



llxs +2 
x 3 = 



3 
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2. f = 1 , t = 1 , h = 1 , g = 1 

3. 

2x 4 -1 

xi = -7 x 4 + 2x 5 + (Ei, i ) 

-3 
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x 4 + 3x 5 



x 2 — 6 + (Ei, 2 ) 

4 

x 5 +2 

x 3 -4 x 5 + (Ei, 3 ) 

3 



4 . 



5 . 

6 . 



7 . 

relation 



2x 4 -1 

xi = -7x 4 + 6y i3 -4 + 

-3 

x 4 + 6yi 3 -4 

x 2 = 6 + (Ei, 2 ) 

4 



2x 4 


+ 


3yn 


-1 = 0 


(Fi,i 


+ 

X 


2x 5 


- 4y 12 


= 0 


( E 1 , 2 




x 5 


- 3yi 3 


+ 2 = 0 


(Fi,3 


Yes . 










Yes : 


|r 35 | = 


1 . Then x 5 = 


3y i3 - 2, the 





(Ti) . Substituting it in the others, we get: 



x 3 



-12yi3+8 + 



3yi3 — 2 + 2 



3 



Remove the equation (Fi, 3 ) . 

Consider f: =2, t: =2; go back to step 3. 



Xi 



2x 4 -1 



(Ei,i) 



(Ei, 3) 



3. 



7 x 4 + 6y 43 -4 + 



-3 



(E 2 ,i) 
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x 4 + 2yi 3 

X2 = Y13+5 + ( E 2 , 2 ) 

4 

X3 = -llyi 3+8 ( E2 , 3 ) 

4. 2x 4 + 3y2i -1 = 0 ( F 2 , 1 ) 

x 4 + 2 yi 3 - 4y22 = 0 (^ 2 , 2 ) 

5 . Yes . 

6 . 

7. Yes. | r 24 I =1. We obtain x 4 = -2y 43 + 4y 2 2, 

relation (T 2 ) . Substituting it in the others we 
get : 

-4yi3+8y22 _ 1 

xi = -2 8y22 + 2 0yi3~4 + (Ea,i) ' 

-3 

X 2 = y 22 + yi3+5 ( h 2 , 2 ) ' 

X3 = - llyi 3+8 (E>, 3) ' 

Remove the equation (F 22 ) • 

Consider f: =3, t: =3 and go back to step 3. 

3. — 2yi3+2y22 _ 1 

xi = -22yi3 - 30y22~4 + (E3 ,i) 

-3 

X 2 = yi3 + y 22 + 5 (E3,2) 

X 3 = — llyia +8 ( E 3 f 3 ) 

4. 2yi3 + 2y22 + 3y3i -1 = 0 (F3 ,i) 
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5 . Yes . 

6 . 

7 . No . 

8. a = 1, and i m =1, ji = 31, j: 2 = 22 and r = 1. 



(A) 


Y 22 — zi - y 3 i. 


relation 


(Hi) . 


(B) 


Substituting 


it in the 


others we get : 








-2yi 3 +2zi-2y 3 i-l 


Xi = 


-22y i3 -30zi + 


30y 3 i-4 + 


(E 3 ,i) 








-3 


x 2 = 


yi3 + z i 


y3i + 5 


(E 3 , 2 ) 


X 3 = 


-Hyi3 


+ 8 


(e 3 , 3 ) 


2y i3 


+ 2zi + y 3 i - 1 


= 0 


(F 3 ,i) 


x 4 


= -2yi 3 + 4zi - 


- 4yi 3 


(T 2 ) 


(C) 


Consider h: = 


2 





10. (B) y 3i = 1 - 2yi 3 - 2zi, relation (Gi) . 

Substituting it in the others we get: 



Xi = 


-40y 13 


— 


92 z 4 


+ 


27 


(E 3 , i ) ' 


X 2 


3yi3 


+ 


3zi 


+ 


4 


(E 3 , 2 ) ' 


X 3 = 


-Hyi3 






+ 


8 


(E 3 , 3 ) ' 


X 4 = 


6yi3 


+ 


12 Zi 


- 


4 


(T 2 ) ' 


y 22 = 


2yi 3 


+ 


3zi 


- 


1 


(Hi) 



Remove the equation (F 3 ,i) . 

Consider g: =2, f: =4 and go back to step 3. 
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Xi = ■ 


-4 0yi 3 - 


92zi 


+ 


27 






(E 4 , i ) 


X 2 = 


3y 43 + 


3zi 


+ 


4 






(E 4 , 2 ) 


X 3 = ■ 


-llyis 




+ 


8 






(e 4 , 3 ) 


No . 


The general 


integer 


solut 


ion of the 


initial 


system is: 














Xi = ■ 


-40ki - 


92k 2 


+ 


27, 


from 


(E 4 .i) 




x 2 = 


3ki + 


3k 2 


+ 


4, 


from 


( E 4 , 2 ) 




x 3 = ■ 


-llki 




+ 


8, 


from 


(E 4/3 ) 




x 4 = 


6ki + 


12k 2 


- 


4, 


from 


(T 2 ) ' ' 




x 5 = 


3ki 




- 


2, 


from 


(Ti) 




where 


ki, k 2 


0 Z. 













Algorithm 4 
Input 

A linear system (1) with not all aij = 0. 

Output 

We decide on the possibility of an integer solution 
of this system. If it is possible, we obtain its general 
integer solution. 

Method 
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1 . 
2 . 



3. 



4 . 



h = 1, v = 1 . 

(A) Divide every equation i by the maximal co- 
divisor of the coefficients of the 
unknowns. If the quotient is not an 
integer for at least one i 0 , then the 
system does not admit integer solutions. 
Stop . 

(B) If there is an inequality in the system, 
then it does not admit integer solutions. 

(C) In case of repetition, retain only one 
equation of that kind. 

(D) Remove all the equations which are 
identities . 

Calculate a = min { I aij | , aij ... 0} and determine 
if j 

the indices io, jo for which this minimum can be 

obtained. (If there are more variants, choose 

one, at random.) 

If a ... 1, go on to step 6. 

If a = 1, then: 

(A) Calculate the value of the variable x-i 

o 

from the equation i 0 ; write this relation 
(V v ) . 

(B) Substitute this relation (where possible) 
in all the equations of the system as well 
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5 . 



6 . 



as in the relations (V v -i) , (H h ) for all v 
and h . 

(C) Remove the equation i 0 from the system. 

(D) Consider v: = v+1. 

Does at least one equation exist in the system? 

(A) If it does not, write the qeneral integer 
solution of the system substituting 

ki, k 2 , . . . for all the variables from the 

right term of each expression representing 
the value of the initial unknowns of the 
system . 

(B) If it does, considering the new data, go 
back to step 2 . 

Write all a± -i, j ... j 0 and b± under the form: 

o o 



O-i j 

0 



= a^ 



cbo 



q v 



+ ri 



with 



r V 



bi 

o 



+ ri 



o o 



with | r± 



o o 



7. Write Xj = - 6 qi j Xj + q ± + t h , relation 

o ■ ■ o o 

>Oo 

(H h ) . Substitute (where possible) this relation 
in all the equations of the system as well as in 
the relations (V v ) , (H h ) for all v and h. 
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8. Consider 

x-i : = t h , h: = h + 1, 

o 

3-i j • — jf j "'jo r 

3-i j • "1 3-i j , bi . L i , 

0 0 — 0 0 0 0 

and go back to step 2 . 

The Correctness of Algorithm 4 

This algorithm extends the algorithm from [4] 

(integer solutions of equations to integer solutions of 
linear systems) . The algorithm was strictly demonstrated 
in our previous article; the present one introduces a new 
cycle--having as cycling variable the number of equations 
of the system--the rest remaining unchanged; hence, the 
correctness of algorithm 4 is obvious. 

Discussion 

1 . The counter variables h and v count the 
relations H and V, respectively, differentiating 
them (to enable the substitutions) ; 

2. Step 2 (A + B) + (C) ) is trivial and is meant to 
simplify the calculus (as algorithm 2); 

3. Substep 5(A) has aesthetic function (as all the 
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algorithms described) . Everything else has been 
proven in the previous chapters (see [4], [5], 

and algorithm 2) . 

Example 4 . Let us use algorithm 4 to calculate the 
integer solution of the following linear system: 



3xi 


- 7x 3 + 6x 4 


= -2 


4xi + 3 x 2 


i 

X 

V) 

+ 


5x 5 = 19 



Solution 

1. h = 1, v = 1 

2 . 

3. a = 3 and i = 1, j = 1 

4. 3 ... 1 . Go on to step 6. 

6. So, 

-7 = 3 ! (-3) + 2 

6 = 3 ! 2 +0 

-2 = 3 ! 0 - 2 

7. xi = 3 x 3 - 2x 4 + ti, relation (Hi) . Substituting 
it in the second equation we get: 

4ti + 3 x 2 + 12x 3 - x 4 - 5 x 5 = 19 

8. Xi : = ti, h: = 2, au.;: = 0, ai 3 : = +2, ai 4 : = 0, 
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an: = +3, b: = -2 
Go back to step 2 . 

2. The equivalent system was written: 

+ 3ti + 3x3 = - 2 

4ti + 3 x 2 + 12x 3 - x 4 - 5 x 5 = 19 

3. a = 1 , i = 2 , j = 4 

4. 1 = 1 

(A) Then: x 4 = 4ti + 3x 2 + 12x 3 - 5x 5 - 19, 

relation (Vi) . 

(B) Substituting it in (Hi), we get: 

Xi = - 7ti - 6 x 2 - 21x 3 + 10x 5 + 38, (Hi) 

(C) Remove the second equation of the system. 

(D) Consider: v: = 2. 

5. Yes. Go back to step 2. 

2. The equation + 3ti + 2x 3 = -2 is left. 

3. a = 2 and i = 1, j = 3 

4. 2 ... 2, go to step 6. 

6. +3 = + 2 ! 2-1 

- 2 = + 2 (- 1 ) + 0 

7. x 3 = - 2ti + t 2 - 1, relation (H 2 ) . 

Substituting it in (Hi)', (Vi) , we get: 
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Xi = 


35ti 


- 6x 2 


- 2 1 1 2 


+ IOX 5 + 


59 




(Hi) ' ' . 


X 4 = 


- 2 Ot i 


+ 3x 2 


+ 12 1 2 ' 


5x 5 - 


31 




(Vi) ' . 


X 3 : = 


ft 

ro 


= 3, 


an : = 


- 1 , a i3 : 


II 

+ 

K> 


bi 


: = 0 


(the 


others 


being 


all = 


0) . Go 


back 


to 


step 2 . 



2. The equation -5i + 2 t 2 = 0 was obtained. 

3. a = 1, and i = 1, j = 1 

4. 1 = 1 

(A) Then, ti = 2t 2 , relation (V 2 ) . 

(B) After substitution, we get: 



Xi = 


4 9t 2 - 


6x 2 + 10x 5 + 


59 


(Hi) " ' ; 


X 4 = 


— 2 8 1 2 + 


3x 2 - 5x 5 - 


31 


(Vi) ' ' ; 


X 3 = 


- 3t 2 


- 


1 


(H 2 ) ' ; 



(C) Remove the first equation from the system. 

(D) v: = 3 

5. No. The general integer solution of the initial 
system is: 

xi = 49ki - 6 k 2 + 10k 3 + 59 
x 2 = k 2 

x 3 = - 3ki - 1 

x 4 = -28ki + 3k 2 - 5k 3 - 31 

x 5 = k 3 
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where (ki, k 2 , k 3 ) 0 Z 3 . 
Stop . 



Algorithm 5 



Input 



A linear system (1) . 



Output 



We decide on the possibility of a integer solution of 
this system. If it is possible, we obtain its general 

integer solution. 



Method 

1. We solve the common system in R n . If it does 
not have solutions in R n , then it does not have 
solutions in Z n either. Stop. 

2. f = 1, v = 1 , h = 1 

3. Write the value of each main variable xi under 
the form: 

(Ef,i) it Xi 0 qijX j + gi + (0 rijX j + ri) / - i, 

j j 

with all qij, qi, rij, r ir -± from Z, so that all 
I r ±j | < | a i|, | r± | < Iml, a i ... 0 (where all Xj ' s 
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of the right term are integer variables: either 

from the secondary variables of the system or 
the new variables introduced with the 
algorithm) . For all i, we write r if j / ~± . 

4. (F f ,i) i : 6 r ± jXj - r if j y f ,i + r± = 0, where 

f 

3 

are auxiliary integer variables. Remove all the 
equations (F f ,i) which are identities. 

5. Does it exist at least one equation (F ffi ) ? If 

it does not, write the general integer solution 

of the system substituting kg, k 2 , . . . for all 

the variables of the right member of each 

expression representing the value of the initial 
unknowns of the system. Stop. 

6. (A) Divide each equation (F ffi ) by the maximal 

co-divisor of the coefficients of their 
unknowns. If the quotient is not an 
integer for at least one io, then the 
system does not admit integer solutions. 
Stop . 

(B) Simplify--as previously ( (A) ) all the 
fractions in the relations (E f ,i)i. 

7. Calculate a = min { | rij | , rij ... 0}, and determine 

i r 3 

the indices io, jo for which this minimum is 




obtained . 



If a ... 

If a = 



1 , go on to step 9. 
1 , then: 




Ill 



(A) Calculate the value of the variable 

o 

from the equation (F f ,i) ; write (V v ) for 
this relation. 

(B) Substitute this relation (where possible) 
in the relations (E f ,i) , (V v +i) , (H h ) for all 
i, v and h. 

(C) Remove the equation (F ffi ) . 

(D) Consider v: = v + 1, f: = f + 1 and go back 
to step 3 . 

9. Write all r± j ... j 0 and r± under the form: 

o o 

ri j = m ! qi i + P j, with | r± j \ < | m | ; 

o o a tr tx 

r± j = ~± ! q± + ri with \r± \ < \ - ± \ . 

oooo o 

10. (A) Write Xj = - 6 q ± j ! Xj + q ± + t h , 

o ■ ■ o o 

D-Oo 

relation (H h ) . 

(B) Substitute this relation (where possible) 
in all the relations (E fri ) , (F f/i ) , (V v ) , 

( H h -i ) . 

(C) Consider h: = h + 1 and go back to step 4. 
The correctness of the algorithm is obvious. It 

consists of the first part of algorithm 3 and the end part 
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of algorithm 4. Then, steps 1-6 and their correctness 
were discussed in the case of algorithm 3. The situation 
is similar with steps 7-10. (After calculating the real 
solution in order to calculate the integer solution, we 
resorted to the procedure from 5 and algorithm 5 was 
obtained.) It means that all these insertions were proven 
previously . 

Example 5 

Using algorithm 5, let us obtain the general integer 
solution of the system: 

3xi + 6x 3 + 2x 4 = 0 

4x 2 - 2x 3 - 7 x 5 = -1 



Solution 

1. Solving in R 5 , we get: 

-6x 3 -2x 4 

Xl = 

3 

2x 3 +7x 5 -1 

x 2 = 

4 



2 . 



= 1, v = 



f 



1, h 



1 
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-2 x 4 

3. (Ei, i ) : xi = 2 x 3 + 

3 

2x 3 +3x 5 -1 

(Ei, 2 ) : x 2 = x 5 + 

4 

4. ( Fi , i ) : -2 x 4 - 3yn = 0 

( F ! , 2 ) : 2x 3 + 3x 5 - 4yi2 -1 = 0 

5 . Yes 

6 . 

7. i = 2 and io = 2, j 0 = 3 

8 . 2 ... 1 

9. 3 = 2 ! 1 + 1 

-4=2! (-2) 

-1 = 2 ! 0 - 1 

10. x 3 = — x 5 + 2yi2 + ti, relation (Hi) . After 
substitution : 

-2x 4 

(Ei, i ) ' : xi = 2 x 5 - 4y i2 - 2ti + 

3 

x 5 +4yi2 + 2ti-l 

(Ei, 2 ) ' : x 2 = x 5 + 

4 

(Fi, 2 ) ' : x 5 + 2ti - 1 = 0 

Consider h: =2 and go back to step 4. 

4 . (Fi,i) ' : -2 x 4 - 3yn = 0 
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(F lr2 ) ' ' : 2ti + x 5 - 1 = 0 

5 . Yes 

6 . 

7. a = 1 and io = 2, j 0 = 5 

(A) x 5 = -2ti + 1, relation (Vi) 

(B) Substituting it, we get: 




-2x 4 



3. 


(E 2 ,i) 


: xi = 


-6ti - 


4yi 2 + 


2 + — 
3 




(e 2 , 2 ) 


: x 2 = 


— 2 1 1 + 


yi 2 + 


1 


4 . 


(f 2 , i ) 


: -2x 4 - 


II 

<N 

>1 

PO 


0 




5 . 


Yes 











6 . 

7. a = 2 and i 0 = 1, j 0 = 4 

8 . 2 ... 1 



9. 



-3 



2 ! ( 1 ) 



1 
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10. (A) x 4 = -y 2 i + t 2 , relation (H 2 ) . 

(B) After substitution, we get: 

(E 2 ,i) ' : xi = -6ti - 4y i2 + 2 + 



2y 2i -2t 2 



(F 2 ,i) ' : -y 2i - 2t 2 = 0 



Consider h: =3 and go back to step 4. 



4 . (F 2 ,i) ' : -y 2 i - 2t 2 = 0 

5 . Yes 

6 . 

7. a = 1 and i 0 = 1, jo = 21 (two, one) . 

(A) y 2i = -2t 2 , relation (v 2 ) . 

(B) After substitution, we get: 



(E 2 ,i) ' ' : xi = -6ti - 4y i2 - 2t 2 + 2 



( H 2 ) ' : x 4 = 3t 2 



(C) Remove the equation (F 2 ,i) . 

(D) Consider v = 3, f = 3 and go back to step 

3. 



3. 


(E 3 ,i) 


: Xi = -6ti - 


4y i2 - 


2 1 2 + 2 




(e 3 , 2 ) 


: x 2 = -2ti + 


yi2 


+ 1 


4 . 
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No . 


The general 


integer 


solution 


of tk 


is : 














Xi = - 


-6ki - 4k 2 - 


2k 3 + 


2 


r 


from 


(E 3 ,i) ; 


x 2 = ■ 


-2ki + k 2 


+ 


1 


r 


from 


(E 3 , 2 ) ; 


x 3 = 


3ki + 2 k 2 


- 


1 


r 


from 


(Hi) ' ; 


x 4 = 




3k 3 , 


from 


(H 2 ) ' ; 




x 5 = ' 


-2ki 


+ 


1 


r 


from 


(V x ) ; 


where 


(ki, k 2 , k 3 ) 


0 Z. 










Stop . 















Note 1 . Algorithms 3, 4 and 5 can be applied in the 
calculation of the integer solution of a linear equation. 

Note 2 . The algorithms, because of their form, are 
easy to introduce in the computer. 

Note 3 . It is up to the reader to decide on the 
algorithm to use. Good luck! 
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